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& angle  of  attack,  degrees 
5 area  of  model,  ft.^ 

\J  free  stream  velocity,  ft. /sec. 

Q suction  quantity,  ft.3/sec* 

A 

free  strean.  dynamic  pressure,  lbs. /ft. 
total  pressure  in  duct,  lbs./ft.^ 

Ho  total  pressure  in  free  stream,  lbs./ft.^ 

Ax  two  dimensional  lift  coefficient 
A*  two  dimensional  drag  coefficient 

£ two  dimensional  pitching  moment  coefficient  about  the  quarter 
”'CU  chord  point 

H - Hi 

A r slot  pressure  coefficient  ■ — ^ * 

-Adi*  two  dimensional  wake  drag  coefficient 

two  dimensional  momentum  drag  coefficient  ■ 
two  dimensional  equivalent  drag  coefficient  ■ /C<% 


■^d  v» 

Ad  c 
Adt 


two  dimensional  total  drag  coefficient  • A*  u>  ■f/Cde. 


Q 


A a suction  quantity  coefficient  » 


Other  engineering  symbols  defined  in  text. 


SUMMARY 


A theoretical  and  experimental  investigation  of  the  primary 
characteristics  of  trailing  edge  suction  is  reported.  The  theo- 
retical study  results  in  a mathematical  method  for  the  prediction 
of  the  performance  of  such  a system  under  limited  conditions.  The 
experimental  study  confirms  the  theory  and  extends  the  investiga- 
tion to  the  determination  of  the  effect  of  suction  slot  geometry 
on  the  lift,  drag,  pitching  moment  and  power  requirement  for  such 
a system.  In  addition,  an  attempt  has  leen  made  to  relate  the 
significant  geometry  of  the  airfoil-slot  comlination  to  the  measured 
lift  increase. 


INTRODUCTION 


One  of  the  operational  requirements  limiting  the  over-all  performance  of 
an  aircraft  is  the  necessity  of  designing  for  reasonably  safe  landing  speeds. 

As  aircraft  maximum  speeds  have  increased  the  speed  range  of  the  airplane  has 
not,  in  general,  kept  pace.  This  has  necessarily  resulted  in  pushing  the  land- 
ing speed  to  the  highest  acccntable  limits  and  perhaps  in  several  instances 
even  beyond  tnese  limits.  Of  course,  this  is  well  recognized  and  designers 
have  made  extensive  use  of  high  lift  devices  to  achieve  even  the  marginal  land- 
ing speeds  that  are  now  accepted.  However,  as  operationax  aircraft  grow  heavier 
and  faster  and  retain  the  same  limits  on  size  (as  in  the  case  of  carrier  lorne 
aircraft)  the  landing  speed  generally  must  increase. 

For  Naval  operations,  the  significant  speed  is  that  of  deck  touch-down 
and  since  this  speed  has  about  reached  its  upper  limit,  a new  approach  for 
safe  carrier  landings  is  desirable.  One  means  would  be  to  build  the  carrier 
larger  and  faster.  The  alternative,  of  course,  is  to  provide  the  airplane 
with  a lower  landing  speed.  Since  compromise  dominates  evary  ahase  of  air- 
plane design  and  operation,  it  is  probable  that  a compromise  will  be  reached 
nere  to  make  the  future  carrier  compatible  with  the  future  naval  aircraft. 

If  present  trends  hold,  wing  loadings  will  continue  to  increase.  This, 
of  course,  suggests  that  a solution  to  the  problem  is  to  find  a means  of  in- 
creasing the  airplane  lift  coefficient  beyond  that  achieved  with  flaps  and/or 
slots.  Several  schemes  designed  to  increase  the  lift  coefficient  of  a profile, 
at  a given  angle  of  attack,  have  been  advanced  over  the  years  and  they  have 
generally  involved  either  blowing  across  some  portion  of  the  wing  or  suction 
at  some  station  on  the  profile.  In  both  systems,  the  boundary  layer  is  af- 
fected, either  energized  with  blowing  or  removed  with  suction.  A combination 
of  blowing  and  suction  has  been  used  and  seems  to  have  some  advantage.  There 
is,  however,  a limit  to  the  amount  of  influence  blowing  or  suction  may  have 
upon  the  boundary  layer.  Beyond  this  limit  lift  continues  to  increase  at  a 
lower  rate,  for  the  case  of  a trailing  edge  suction  slot,  seemingly  due  to 
the  direct  influence  of  the  sink  upon  the  adjacent  free  streamlines.  For  this 
reason,  the  term  circulation  control  or,  more  recently,  super  circulation  has 
been  applied  to  these  devices.  A lift  increase  caused  by  the  removal  of  some 
portion  of  the  boundary  layer  is  equivalent  to  an  increased  circulation.  There- 
fore, strictly  speaking,  it  seems  improbable  that  the  two  types  of  lift  increase 
can  be  separated  and  the  system  is  finally  one  of  increasing  circulation  by  one 
means  or  another.  Thus  small  amounts  of  suction  remove  the  boundary  layer,  or 
some  portion  of  the  boundary  layer,  thereby  permitting  the  bordering  streamlines 
to  follow  more  closely  the  pattern  of  potential  flow  resulting  in  increased  cir- 
culation. Larger  quantities  of  suction  seem  to  influence  the  adjacent  stream- 
lines directly,  inducing  a greater  downward  component  to  the  flow  which  again 
provides  greater  lift  and,  of  course,  more  circulation. 

Princeton  University  became  interested  in  this  phase  of  low  speed  flight 
research  with  a suggestion  from  the  Office  of  Naval  Research,  Air  Branch,  that 


the  validity  of  the  early  Regenscheit  work,  in  the  field  of  boundary  layer 
control  by  trailing  edge  suction  slots,  be  investigated.  This  was  formal- 
ized by  Contract  Number  N6-onr-270l6  in  March  of  1951  between  the  Office  of 
Naval  Research,  Air  Branch  and  Princeton  University,  to  conduct  a general 
investigation  of  trailing  edge  suction. 

The  program  of  investigation  subsequently  setup  was  both  extensive  regard- 
ing slot  shape  and  location  and  intensive  with  respect  to  the  thousands  of 
individual  measurements  of  lift,  drag,  pitching  moment,  suction  quantity  and 
power  requirement  of  the  various  configurations  tested.  Since  the  effects  of 
the  several  variables  of  a trailing  edge  suction  slot  were  largely  unknown 
quantities,  it  was  decided  to  limit  the  stucty  to  the  airfoil  without  a flap. 
Future  experiments  are  to  be  made,  however,  which  will  include  a flap  study. 

This  report  deals  with  the  problem  of  trailing  edge  suction  both  theo- 
retically and  experimentally.  The  principal  aim  has  been  to  understand  and 
evaluate  the  mechanics  of  the  suction  slot  geometry.  Since  any  system  of 
circulation  control  must  be  evaluated  with  respect  to  a given  design,  the 
results  reported  here  should  be  helpful  in  such  an  evaluation. 

XI  Reference  1. 
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I.  THEORY  OF  TRAILING  EDGE  SUCTION 


In  order  to  understand  the  basic  mechanics  of  increased  circulation 
due  to  trailing  edge  suction,  it  is  first  necessary  to  examine  the  pos- 
sible flow  patterns.  When  these  flows  are  understood,  it  becomes  possible 
to  arrive  theoretically  at  an  expression  relating  A sCjl  * and  a 

suction  slot  pai'ameter. 

A.  Determination  of  the  Type  of  Flowt 

A study  of  the  flow  pattern  necessarily  requires  that  considera- 
tion be  given  to  the  location  of  the  stagnation  points.  In  general,  it 
may  be  stated  that  a stagnation  point  leaves  the  lower  lip  of  the  trailing 
edge  for  the  free  stream  if  the  suction  coefficient  exceeds  a certain  value. 
This  can  be,  and  has  been,  visualized  in  a smoke  tunnel  study.  The  lift  and 
the  suction  quantity  at  the  instant  the  stagnation  point  leaves  the  profile 
can  be  determined  theoretically.  This  point  may  be  considered  the  combina- 
tion of  two  stagnation  points,  one  of  them  being  located  on  the  trailing 
edge  satisfying  the  Kutta-Joukouski  condition,  and  the  other  located  between 
the  sink  and  the  trailing  edge  (See  Figure  1-a  and  1-b).  The  stagnation 
point  T moves  aft  with  increasing  suction  until  it  reaches  the  trailing 
edge  where  it  is  combined  with  7 It  • Any  increase  in  suction  coefficient 
forces  this  combined  stagnation  point  to  the  free  stream  (see  Figure  1-c 
and  1-d). 

Figures  2-a,  2-b,  and  2-c  illustrate  flow  conditions  that  are 
theoretically  possible  and  predictable  under  certain  circumstances.  In 
these  T continues  to  move  around  the  trailing  edge  to  the  lower  surface 
of  the  profile.  The  case  is  shown  (Figure  2-b)  where  T meets  the  stag- 
nation point  T,  moving  aft  on  the  under  surface.  An  increase  of  suction 
coefficient  at  this  point  causes  this  combined  stagnation  point  to  jump  to 
the  free  stream  as  shown  in  Figure  2-c.  The  motion  of  the  stagnation  points 
can  be  studied  analytically  in  the  following  way.  The  function 

H * fft) 


yields  a one  to  one  conformal  representation  of  the  area  /«?/>/  into 
the  outside  of  the  profile  (in  the  z-plane)  in  such  a way  that 

so  that  I corresponds  to  the  trailing  edge  of  the  profile.  This  flow, 
having  a sink  at  the  point  corresponding  to  , can  be  represented 

by  the  complex  potential  function  ' 


x = r + *r)+  & ~ € 


a) 


For  a stagnation  point  i is  valid 


For  the  stagnation  point  on  the  trailing  edge  $ « / 

^ e‘«>  - # 7^>  = o 

Subtracting  (3)  from  (2)  yields: 


Dividing  thru  by  J - / results  in 

^o€to<-^r  ~ “ 

Multiplication  by 

(u/s  e‘*-  J r* + 

It  follows  from  (3)  that 
Substituting  in  (U)  yields 

[^oeoM(i-et6)  + $ jf^zs — ^ * *- c 

' C * 77^  T^7*'  ) S ~ ^ =0 

This  may  be  re-written  as 

p-»  A 
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r~ 


In  addition  to  the  stagnation  point  /•»/  there  are  two  other  stag- 
nation points  as  solutions  of  equation  (E>)»  By  denoting 


the  solutions  become 


3=';c”*')(-<s±VTr^) 


(7) 


For  fixed  values  of  »< 
the  amount  of  suction  Q.  This 


• If  Isle  I,  then 


, & ,0/5  the  values  of«Sare  variable  with 
variation  is  such  that  J increases  with  Q 


/cs  = i 


(8) 


which  means  that  both  stagnation  points  are  located  on  the  circle  )^/=|  . 

To  these  correspond,  due  to  , two  stagnation  points  on  the 

profile  surface.  If,  however,  IS/*/  then 

or  . (9> 

One  of  these  points  is  outside  the  circle  (see  Figure  3)  I Si  “ ty  on  the 
radius  which  is  inclined  at  the  angle  c*  >+  /3/e  from  the  negative  axis  , 
the  distance  being  S+  V SA- 1 from  O . The  other  uoint  is  situated 

inside  the  circle  on ‘the  same  radius, the  distance,  in  this  case,  being 

S-V5^7  from  0.  Only  the  first  of  these  is  of  interest  since  it  corresponds 
to  a stagnation  point  outside  the  profile  in  the  z-plane,  With  increasing 
suction  (Q)  this  point  moves  toward  infinity  in  the  free  stream  and  the  con- 
dition for  its  existance  is 

lsl>/  (10) 

Obviously,  the  determination  of  which  of  these  flow  patterns  can  be 
expected  is  of  basic  importance  to  the  under standing  of  trailing  edge  suc- 
tion. A smoke  tunnel  experiment  (See  Figure  ii)  clearly  proved  that  the  free 
stream  stagnation  point  of  Figure  1-d  is  the  solution  with  physical  signifi- 
cance. With  a stagnation  point  of  this  type, flow  around  the  lower  trailing 
edge  occurs  which  greatly  reduces  the  effect  of  trailing  edge  suction*  Hence, 
unless  measures  are  taken  to  prevent  this  motion  of  the  stagnation  points,  the 
value  of  suction  coe  ficient  which  causes  the  two  s tagnation  points  to  combine 


*r~ 


u. 


at  the  trailing  edge  can  be  considered  the  practical  upper  limit  for  this 
means  of  increasing  lift. 

B.  Lift  Increase  due  to  Trailing  Edge  Suction: 

The  theoretical  developments  presented  in  this  section  are  based 
upon  a normal  wing  profile  with  a small  suction  slot  located  on  the  upper 
surface.  It  will  be  assumed  that  this  slot  has  parallel  straight  walls  and 
that  suction  is  produced  by  a sink  6 situated  within  the  slot.  To  satis- 
fy continuity  requirements,  a source  of  equal  strength  must  be  assumed  to 
be  located  at  infinity.  It  is  further  assumed  that  the  profile  is  at  an 
angle  of  attack  <x  in  a uniform  flow  field  having  the  velocity  ^at  in- 
finity. Figure  5 indicates  the  resulting  flow  field  for  a moderate  suction 
quantity.  For  this  suction  quantity  a stagnation  point  P appears  on  the 
upper  surface  of  the  profile  between  the  suction  slot  and  the  trailing  edge. 
(iMs  is  also  illustrated  by  Figure  6).  The  flow  leaves  the  trailing  edge 
smoothly,  indicating  the  presence  of  a stagnation  point  attached  to  the  pro- 
file at  this  point  and  hence  the  Kutta-Joukowski  condition  is  satisfied  as 
in  the  case  of  no  suction.  A third  stagnation  point  is,  of  course,  situated 
near  the  nose  of  the  profile. 


The  two-dimensional  potential  flow  of  an  incompressible  fluid  satis- 
fying these  conditions  can  be  determined  in  the  classical  manner.  To  do  this, 
the  outside  of  the  profile  in  the  t ~ * * *•  $ plane  is  transformed 
conf  ormally  to  the  outside  of  a circle  of  unit  radius  in  the  J * Jf  ■+  <■ 
plane  in  such  a way  that  the  infinite  points  of  the  two  planes  correspond 
and  that  the  size  and  position  of  the  profile  in  the  9 plane  is  such  that 
d#S  *1  for  and  JT»  / corresponds  to  the  trailing  edge  of  the 

profile.  (See  Figure  5).  Under  these  conditions,  the  transformation  is 
determined  uniquely  and  the  region  outside  of  the  circle  I Si- I has  an 
expansion  of  the  foim 


J+-f  • ■ 


(ID 


Under  the  restrictions  of  the  transformation,  the  x-axis  corresponding  to 
O(»o  is  at  the  same  time  the  axis  of  zero  lift  without  suction.  The  length 
of  the  profile  is  determined  by  the  transformation  and  has  a value  approxi- 
mately equal  to  four. 

* 0 

As  a result  of  the  transformation,  a well  defined  point,  o * “ 
where  & is  the  arc  of  the  circle  between  the  points  corresponding  to  the 
trailing  edge  and  the  slot,  corresponds  to  the  slot  exit  on  the  profile. 

Kore  rigorously,  the  point  determined  by  corresponds  to  the  sink  within 
the  slot,  and  the  slot  walls  correspond  to  two  small  arcs  adjacent  to  fi  if 
the  inside  of  the  slot  is  included  in  the  transformation.  If,  however,  the 
slot  is  considered  infinitesimally  small,  fi  defines  a point  which  corresponds 
to  the  slot  including  the  sink  as  a limit. 


The  flow  under  discussion  is  determined  by  the  complex  potential  function 

X = “A  (Se1'**  7^**)  + X + 


(i?) 


5 


The  first  tern  represents  the  parallel  flow  at  an  ancle  of  attack 
OC  around  the  circle,  the  second  term  a circulation  (vortex)  around  the 
circle  with  strength  T1  . The  third  term  results  from  a source  with  the 
strength  Q at  Js^,  a source  ,of  equal  strength  at  3=  and  sink  with 
the  strength  -a<?  at  3 = so  that  the  strength  of  this  sink  with 

respect  to  the  flow  outside  the  circle  is  equal  to  -Q,  where  Q is  measured 
in  units  of  volume  per  unit  span  and  time.  The  source  at  3 s ° inside 
the  circle  is  added  in  order  to  obtain  the  contour  of  the  circle  as  a stream- 
line. The  flow  represented  by  is  transformed  by  i*  ■pa) 

to  a flow  with  corresponding  singularities  in  the  plan  of  the  profile. 

The  complex  velocity  function  may  be  expressed  as 

dx  _ 


f = ure  1 


(13) 

where  'fT3)=  ,u r is  the  velocity  and  Q is  the  direction  of 

the  flow  at  a point  in  the  X plane.  Thus,  for  a stagnation  point  in  the 
* plane  s O 

Hence^fnr  a point  JT^C  c * to  correspond  to  a stagnation  point  on 

the  profile^  ■H-  must  equal  zero  at  the  point^with  the  further  condition 
that  -f'(£)=  if  is  not  zero  or  is  at  least  of  smaller  order  than  • 

Differentiation  of  equation  (I2)yields: 

7*-)+  f + & (j-  ~ 

ir 


~df=  C'o  I 


For  a stagnation  point  = on  the  circle  1*1*1 

sary  that  * — Q , Therefore; 


it  is  neces- 


usbfe  - 

Multiplying  by 


3 

‘.(rni-Jg)' 


O+jFr?  *'*  + & (***-  )=  ° 


y 


or 


1 -C  (y— 

1.  Pc 

■e  J 

+2.7T 

*) 

l.Ct 

-«•  ; 

'^77- 

‘C«-8)  “« 

— C 

; fo<-  a 

cY(B-0 


= 2 i (*  -£r) 


now 


ii 


6 


and 


Therefore 


g ± e-  ^ 
'W~rrW 


— - L COT 

SL 


r=  -4  TTOJ-o  •s/'/v<V-«)  + Q£»~r  ^C 


(lfc) 


where  the  angle  ^ determines  the  position  of  any  stagnation  point.  Since 

the  point  fs/  (\C.  V—  o')  is  supposed  to  be  a stagnation  point  in  accord- 

ance with  the  Kutta-Joukowski  condition,  the  relationship  simplifies  to 


r*=4  1TUZ  $/*«  + QCcrr%L  ( 1 $) 

Thus  the  circulation  is  determined  by  the  velocity  usp  , the  angle  of  attack 
o(  , the  position  (9  and  the  strength  Q of  the  sink. 

The  suction  coefficient  is  generally  defined  as 


sCq  — 


-A- 

u /is 


(16) 


where  3 is  the  wing  area,  UJb  the  free  stream  velocity  and  Q the  volume  of 
suction  per  unit  time.  In  the  preceeding  discussion,  a unit  span  has  been 
assumed,  thus  if  G is  the  profile  chord  length,  /C  Q may  he  defined  as 


/•  , — & — 

Q — C 


(16a) 


where  C is  approximately  equal  to  four.  If  T is  now  eliminated  from 
equation  and  16,  one  obtains  the  coefficient 

-4-Tr *•'*  fa- Sl - 
X.G)  — C co-rApc — cot  Sr 

which  simplifies  to 


XTo=  -^‘cosfa  — (17) 

From  this  relation  one  can  determine  how  the  stagnation  point  between  the 
suction  slot  and  the  trailing  edge,  defined  by  Jf^moves  with  changing.,^  $ . 

As  previously  explained,  a case  of  special  interest  occurs  when/Cq 
is  increased  sufficiently  to  cause  this  stagnation  point  to  reach  the  trail- 
ing edge.  In  this  case  reaches  the  limit  zero  and  one  obtains  the  value 

;€<**  = -i^C0S«  ■ (18) 

where  is  the  highest  suction  coefficient  for  which  trailing  edge 
suction  is  practical. 
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The  components ? X and  y of  the  force  acting  on  the  outside  of  the 
profile  due  to  the  pressures  produced  by  the  flow^can  be  det' mined  from 
Blasius*  equation? 

''  Y+  * x = - « £ (&)*<>•  (19) 

where  ^ is  the  density  of  the  fluid.  In  this  particular  case,  the  curve 
of  integration  C which  must  be  considered,  is  the  surface  contour  of  the 
profile  up  to  the  slot,  the  parallel  sides  of  the  slot  and  a perpendicular 
connection  between  the  sides  of  the  slot.  When  transformed  into  the  $ -plane 
equation  19  becomes 


ds  <20) 

The  curve  of  integration  K now  consists  of  the  circle  l$l-  I up  to  a 
small  half  circle  excluding  the  point  S=  from  the  outside  region 

and  of  this  small  half  circle  itself  which  corresponds  to  the  connection 
of  the  slot  walls  . This  curve  R can  b e deformed  into  a large  circle 
around  because  there  are  no  singularities  of  the  integrand  outside 

except  those  at  infinity. 

Differentiation  and  simplification  of  equation  12  yields: 


d£  - (o-L«  £ll)  + SLf  rv  jL  _ L 

^ + 277  1 7r  % - 

4 

Expanding  this  expression  into  a power  series  of  using  the  expansion 

rt1* : fO  * f" + ~rt~  *■  ■ ■) 


results  in 


dff  = J ■ 


and 


a f-)we-*‘v  - ax*-  ■ ■ 

The  function  has  the  expansion 


8. 
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Equation  20  yields  then 

['*$*■  *■■ 
or 

r+tf*  ■ ■ ] ds 

or  '• 

V+4-  x=  <°«Ji(r+Q<)e-~‘'*  (21) 

The  components  of  this  force  parallel  and  perpendicular  to  the  free 
stream  direction  are  the  drag,  d , and  the  lift, & j therefore, 

Bfy'+tXjs4* 

= furo(r+Qc) 


Thus 

, d=^c^,o  (22) 

This  result  has  been  derived  by  C.  B.  Smith^.  Substituting  the  value  of 
T*  from  equation  15  into  equation  22,  one  obtains 

•4=  fur0  (4-7Tcor0  s/^a  + q cot  P/l)  (23) 

By  definition  the  lift  coefficient  is 

~ /%  ^ 
Substituting  equations  23  and  16-a  in  equation  2h  yields 


sCjl  = + £/Cg(  cot  fyz  W) 

An  examination  of  equation  (2$)  reveals  that  the  lift  coefficient  con- 
sists of  ^“sinot,as  predicted  for  the  zero  suction  case  by  classical 
aerodynamic^  plus  a term 

A/Cjl  = Z/C^  Ccrrfi/z  (26) 

due  to  suction.  When  the  stagnation  point,  ~P  , reaches  the  trailing  edge, 
the  expression  for  the  lift  increase  becomes 


1.  See  Reference  2. 


1 p t 
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ACj*  - cTT  C (27 ) 

which  represents  the  upper  practical  limit  for  trailing  edge  suction. 

C.  The  Determination  of  (3  in  the  Case  of  Trailing  Edge  Suction. 

If  the  suction  slot  is  small  and  enters  the  surface  of  a profile 
with  parallel  walls,  the  value  of  |B  can  Le  found  by  the  conformed  mapping 
of  the  profile  into  a circle  as  previously  described.  If,  however,  as  is 
usually  the  case  of  trailing  edge  suction,  the  end  of  the  profile  has  been 
cut  off  so  that  the  edfees  of  the  upper  and  lower  surfaces  form  the  slot 
entrance,  then  the  conformal  mapping  of  the  original  profile  will  not  lead 
to  a correct  determination  of  p . This  arrangement  has  now  to  be  treated 
as  an  entirely  new  profile  like  that  of  Figure  7,  the  location  of  the  real 
slot  being  indicated  by  an  arrow. 

For  this  condition  the  following  procedure  is  proposed.  The  edge 
of  the  longer  lip  of  the  trailing  edge  (usually  the  lower  lip)  is  connected 
by  a straight  line  to  the  mid  point  of  the  radius  connecting  the  center  of 
curvature  of  the  leading  edge  to  the  leading  edge.  This  straight  line  is 
then  considered  a slot  and  the  region  surrounding  it  is  transformed  con- 
formally to  the  region  outside  a circle.  If  the  straight  line  connects  the 
points  2*  — 7?  and  M + t"R  in  the  plane  of  the  profile,  the  transformation 
to  be  applied  is  the  inverse  transformation  of 

zm  y*(5+ -f^)  (28) 

(the  Joukowski  transformation)  or 

J = (29) 

This  transformation  is  frequently  applied  in  order  to  transform  a profile 
into  a nearly  circular  curve.  The  numerical  compulation  of  f*  { ^ as  a 

function  of  2 = X-fc/  can  be  performed  by  expressing  J in  terms  of 
polar  coordinates 

t <p 

or  j=re 

and  substituting  this  in  equation  28  yielding 

* » 'A  (Ve^+  ^rCc*) 

By  separating  the  real  and  imaginary  parts 

x s '/e.(>r+¥)cos4>  } (30) 

Y ~ J 

Substituting  in  equation  30  for  V* 

r 


where  S is  a new  variable, results  in 

* =■  7? 5 cos  ft 


Y-  ~Rs/A/h  6 s/A/<f> 


and 

*± , X— = / 

K*Cosh*G  "** 

Replacing  in  equation  31  by  {CoSh^G 

algebraic  equation  of  second  order  in  (LOSh*  O . 
equation  yields 


(31) 

~ I ).  one  obtains  an 
The  solution  of  this 


CoSA~&  - £ +l)+  \jfa  (— ftV  * 


(32) 


Knowing  Cos/ , O,  the  relation 

• C OS  (ft  = -ftCoSAO 


(33) 


may  be  applied  and 


y=  ~R  ee 


(3« 


then 


f -tiCosft  (35) 

Figure  8 shows  the  transformed  profile  points  of  figure  7 where 
the  trailing  edge  of  a NAO^  230l£  airfoil  has  been  cut  off  10£  on  the  upper 
side  and  $%  on  the  lower  side.  The  transformed  profile  is  now  approximated 
by  two  tangential  circular  arcs,  one  ending  at  the  upper  edge  A and  the 
other  ending  at  the  point  S-  "R  . Both  arcs  are  tangent  at  a point  C some- 
where near  the  point  corresponding  to  the  profile  leading  edge.  A third 
circle^  orthogonal  to  these  two;  passing  through  A and  C is  now  considered. 
This  circle  (in  this  example  a straight  line)  intersects  the  upper  arc 
RC  at  the  point  B.  The  arc  AB  now  corresponds  to  a curve  in  the  z-plane 
which  closes  the  profile  by  connecting  the  edge  of  the  upper  lip  with  a 
point  on  the  line  between  -R  and  + R forming  the  upper  boundary  of  the  lower 
lip.  This  connection  is  perpendicular  to  both  lips  and  can,  as  an  approxi- 
mation, be  replaced  by  a small  circular  arc.  On  this  connection  the  slot 
entrance  S (sink)  with  parallel  sides  is  assumed  as  indicated  by  the  arrow 
in  Figure  8.  This  configuration  corresponds  very  closely  to  the  normal 
geometry  of  a profile  with  trailing  edge  suction  . 

A transformation  by  reciprocal  radius  with  C as  the  center  is  now  used 
to  map  the  outside  of  the  area  bounded  by  the  three  arcs  into  a half  plane 


with  a rectangular  step.  (See  Figure  9)  This  stepped  half  plane  can  be 
conformally  transformed  by  use  of  the  Schwartz-Christoffel  transformation 
into  a half  plane.  If  + is  the  complex  coordinate  in  the  half 

plane  as  indicated  by  Figure  9,  the  transformation  is 


t?o  and  t « — / are  the  points  corresponding  to  the  corners  of  thestep. 
Integration  of  this  expression  yields 


5=  +j£y>(/T 

From  this,  one  secs  that  3 = 0 corresponds  to 
T.  -/. 

In  order  to  compute  this  transformation,  it  is 
a new  complex  variable  T defined  by  the  relation 


(36) 

dvtd  S-JF c to 

convenient  to  introduce 


fE  + fT+t  ~ T 


(37) 


Sot 


ft  * / V/v-r  */<(T-h+) 


and 


(38) 

s*y+  (rM-  ■£)  ^ r 

By  letting 

T=  ere tC 

(39) 

one  obtains  the  expressions 

S , * -atycos  e€+J2>iO- 

(l*Oa) 

S*.-  '/4-(o-*  + rX)Si  * 

(liOb) 

and 

t,  = t ^2)C0S  Z€ 

(lOa) 

*a= 

(lab) 
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The  values  of  S corresponding  to  real  values  of  t 2 O can  be  found 
directly  from  equation  36  and  are  plotted  in  Figure  10.  In  order  to  ob- 
tain the  values  of  <r~  which  correspond  to  points  on  the  straight  line 
3,  so,  equation  b Oa  must  be  used  to  compute  € as  a function  of  a~  . 

The  other  three  expressions,  equations  bOb  and  Ula  and  b yield  the  values 
SAi  trt  and  as  functions  of  the  same  parameter,  . Figure  11 

presents  value0  of  r,  and  plotted  as  functions  of  J*  • This  figure, 
together  with  Figure  10, contains  sufficient  information  for  the  determina- 
tion of  & . 

The  t-half  plane  is  now  transformed  into  the  region  outside  of  a 
circle  in  such  a way  that  infinity  of  the  z-plane  corresponds  to  the  in- 
finite point  of  the  now  plane.  This  can  be  done  by  a transformation  by 
reciprocal  radius  with  C'  as  a center,  where  C'  is  the  point  corresponding 
to  C after  the  step  transformation  (Figure  8).  In  this  manner,  the  real 
*tr  axis  is  transformed  into  a circle  K . 

The  two  inversions  and  the  s tep  transformation  have  to  be  constructed 
or  computed  for  only  the  original  points  R and  S (corresponding  to  the  sink). 
If  R ' and  S / are  the  images  on  the  circle  K with  the  center  K,then 

For  the  example  of  Figure  7 where  the  end  of  the  NACA  23015  airfoil 
has  been  cut  off  from  90%  chord  on  the  upper  surface  to  95%  on  the  lower, 
the  above  method  yields 

A/Cx  - Z, C q Cor  /%  r /A.SsCq. 

In  Figure  13,  this  theoretical  result  is  compared  with  experimental 

Pressure  distribution  measurements  for  •<»©*.  The  actual  t railing  edge 
b used  in  these  measurements  is  shown  in  Figure  12.  If  trailing  edge  la 
of  this  figure,  where  the  lower  lip  reaches  the  trailing  edge  and  the  up- 
per lip  is  maintained  at  90%  chord,  is  used,  the  theory  yields 

A/Cjx  = 6.0  SsC  $ 

A balance  measurement  foro(*c0is  compared  for  this  case  with  the  theo- 
retical result  in  Figure  lb.  In  both  cases  the  maximum  AA  values  at- 
tained are  considerably  lower  than  the  corresponding  A values,  from 
which  it  may  be  concluded  that  the  s tagnation  point  P did  not  reach  the 
trailing  edge  in  these  tests. 

D.  Computation  of  the  Pressure  Distribution! 

The  methods  described  in  the  previous  section  may  also  be  used 
to  compute  the  pressure  distribution.  In  this  case,  all  the  points  of 
the  profile  transformed  by  the  inverse  Joukowski  transformation,  shown  in 
Figure  8,  must  be  further  transformed  by  the  reciprocal  radius  and  step 
transformations.  These  transformations  yield  a set  of  points  in  the 
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neighborhood  of  the  circle  K.  By  any  of  the  well  known  transformations, 
these  points  may  be  transformed  into  points  on  the  circle  K itself. 

A simpler  and  as  accurate  a method  is  to  use  the  conformal  trans- 
formation of  the  original  profile  and  to  locate  a suction  slot  at  the 
point  corresponding  to  S ~ as  found  by  the  preceding  method.  It 

is  advisable  to  use  this  method  particularly  if  the  conformal  transfor- 
mation of  the  original  profile  is  known  as  in  the  case  of  many  of  the  NACA 
profiles . 


In  order  to  find  the  velocity  of  the  flow,  the  expression 

"•m 

may  be  computed  by  differentiating  and  simplifying  equation  12,  yielding 

ck  - Cat  (*U  -Of  ' <3  I 

dS  -or0(e  - zit  ~ tt  Jr 


The  values  of  the  points  on  the  profile  correspond  to 

JT  = e.1* 


Thus 


it  - 

multiplied  by  e ‘ ^ results  in 

. i<P 


>)  cf*-#))  Q+ri. — 12. — — 

— c htt  t 

- a/-.  i (•«.«)  zx  4 _a/;_ 


Now 


' ~ e‘^-  e4*  £*•  fl c ~ 

Therefore 


i 0**2- 
e*v_  - V4  *• 


-2LUr°  1 %■  <-Coe  tjP 

Replacing  P by  Q using  equation  l£  yields 

= w(?b- + -$r(cot  ^ ^ f Co  r 

By  introducing  Cq  (from  equation  16a)  in  place  of  Q in  the  above  expression 
and  by  letting  C<4  as  an  approximation,  this  equation  yields 

|g^|  = U/S  |£(V/V(0-<X)f  £//V«)f  4-)| 


Ill 


T 


Now 


urz\±)L 

\d* 


AJL 


d$ 


Combining  equations  h2  and  h3  results  in 

- + *£*(c-ot  4-J*-  + cot 


ur  ~ 


(U3) 

m 


In  this  equation  <p  is  the  amplitude  of  the  point  Js-e4-^  corresponding 
to  the  profile  point  under  consideration.  It  should  be  noted  that  the 
equation  has  been  developed  on  the  asiumption  that  the  chord  C is  equal  to 
. Using  this  assumption,  it  is  possible  to  express  the  ordinate  as 

X = 2.  COS  <p 

with  only  a small  error.  If  the  velocity  distribution  of  the  profile 

is  known  for  the  case  without  suction  {/Co.  ■*  O ) for  a given  angle  of  attack 
o<  , the  value  of  \f'( $)j  can  be  found  for  any  point  X~2.cos  <p  on  the  pro- 
file by  use  of  equation  hiu  This  equation  can  then  be  used  to  obtain  the 
velocity  distribution  for  any  q (3  and  c*  . The  pressure  dis- 

tribution is  then  obtained  from  the  equatior 

In  this  manner,  the  computed  pressure  distributions  of  Figures  1$  and 
16  have  been  obtained  for  comparison  with  the  experimental  measurements. 


II.  THEORY  OF  SUCTION  SLOT  DESIGN. 


The  magnitude  of  the  lift  produced  by  suction  through  a slot  is  mainly 
dependent  on  the  position  of  the  sink  which  replaces  the  slot.  However,  the 
economy  of  the  lift  producing  device,  assuming  a fixed  position  of  the  sink, 
will  depend  on  the  shape  of  the  suction  slot.  In  order  to  produce  a laminar 
flow  within  and  near  the  slot,  it  is  necessary  to  determine  the  correct  slot 
shape  30  that  there  is  no  decrease  in  flow  velocity  along  the  slot  walls  in 
the  direction  of  the  flow.  The  hodograph  method  is  used  to  determine  this 
shape  for  the  case  of  a trailing  edge  suction  slot. 

A.  Analysis  of  the  Problem! 

The  trailing  edge  of  the  original  profile  is  assumed  to  be  formed 
by  two  planes  including  the  angle  (Figure  17).  This  trailing  edge  is 
cut  off  so  that  the  lower  surface  has  the  end  point  0=C  and  the  upper 
surface  the  end  point  A.  The  trace  of  the  slot  in  the  x-y  plane  is  formed 
by  the  straight  line  ending  at  o=C  as  the  lower  wall  and  an  undetermined 
wall  connecting  point  A with  point  B,  the  latter  being  the  theoretical  in- 
finite point  of  the  lower  wall.  The  flow  entering  the  slot  will  have  a 
stagnation  point  S on  the  lower  wall  and  will  leave  the  edge  o = C smoothly 
according  to  the  Kutta-Joukowski  condition  which  is  valid  in  this  case. 


The  wall  AB  has  to  le  designed  in  such  a way  that  the  flow  has  increasing 
velocities  from  K to  A,  from  S to  B and  from  S to  o*C.  and  a constant  velo- 
city along  the  wall  AB.  These  conditions  are  imposed  in  order  to  prevent 
flow  separation  s.'nce  it  is  known  that  avoiding  an  adverse  pressure  grad- 
ient prevents  separation  on  a flat  plate.  It  is  by  no  means  clear,  however, 
that  the  avoidance  of  an  adverse  pressure  gradient  has  the  same  effect  in 
regions  of  high  curvature.  In  Figure  17  the  flow  directions  are  indicated 
by  arrows  and  numbers.  The  slot  width  is  designated  by  6 and  K denotes  the 
point  on  the  upper  surface  where  the  flow  has  the  lowest  velocity.  The 
sink  is  located  at  B.  The  function 

= JL>*.  -g-g  = US-+L& 

represents  the  logarithmic  hodogr^ph,  ?c  being  the  complex  potential  func- 
tion, w the  flow  velocity  and  ^ the  direction  of  the  flow  for  the  point 
xt  — x + i 3 . The  flow  in  the  z-plane  corresponds  then  to  the  flow 

represented  by  Figure  18  in  the  t plane.  The  notations  for  corresponding 
points  and  flow  directions  are  the  same  as  in  the  z-plane.  The  constant 
velocity  along  the  wall  AB  is  denoted  by  , the  finite  velocity  at  C 
by  ur%  and  the  velocity  minimum  at  K by 

A one  to  one  conformal  tr^nsf ormation  of  the  flow  area  in  the  z-plane 
into  the  upper  half  of  the  3 - £ + *■  >1  plane  is  represented  by  Figure  19. 
Again,  the  corresponding  ooints  and  flow  directions  are  denoted  by  the  same 
symbols.  The  transformation  ic  accomplished  in  such  a way  that  3 — °ri 
corresponds  to  the  infinity  outside  the  trailing  edge  in  the  z-plane  and 
that  the  points  S-  -/  and  correspond  respectively  to  the  points 

A and  E.  The  abscissas  Ji  of  the  real  points  corresponding  to  K,S,  and  C 
are  denoted  by  m,  s,  c.  The  flow  in  the  JT  -plane  las  a source  at  infinity', 
a sink  at  S 3 / and  stagnation  points  of  first  order  at  3 ■»  ^ and 
S-C  . 

B . Transformation  of  the  3 -Plane  Into  the  t-Plane: 

The  flow  area  in  the  t-plane  is  a polygon.  The  conformal  trans- 
formation of  the  upper  J half  plane  into  this  polygon  can  be cfetermined, 
therefore,  by  Schwartz-Christoffel’s  equation.  This  yields 

*•-  * c.,  ^ 

where  C0  and  C,  are  constants  to  le  determined. 


(U6) 

Starting  with  a large 
t is  augmented  by 

Co  7Tc  = 7TL 


~£~  — (2o  j2vi 


The  integral  con  be  evaluated  ai}d  gives 

(Ss- 

The  logarithm  is  real  for  real  values  of 
positive  value  of  3 and  massing  over  S the  function 


+ C, 


16 


r~ 


therefore 


c 

Passing  with  $ over  3 

l / &*■- 1 
S-/*7 


_ j {lEZL 

o - s-/y> 


= + oo  the  function  t is  augmented  by 
7Tl  -^c. 


therefore 


a.  , £EZ 

7T  5 -/>7 


(U7) 


(U8) 


An  inspection  of  equations  U7  and  iih  shows  that 

Ca-~§T  (h9) 

Thus  C*  is  determined  by  @ . Assuming  a given  value  for  S,  it  follows  from 
equation  U8  that 

/n-  s - !/•$*■-  / (5o) 

So  the  position  of  the  velocity  maximum  is  determined  too. 

Approaching  S-+I  along  the  real  S axis  from  the  right  side,  we 

find 


* - 7TL  C,  ur,  y-  7 ri 


Multiples  of  2 7Tt.  to  be  added  to  t eventually  are  without  importance 
because  =■  . Thus 


61/  ~ urt 

and  the  result  is 


(51) 


t = Jtv 1 


Ts  -/  J 


(52) 


valid  for  /*/  =*  /• 

For  l$l=/  equation  U5  is  replaced  by 

t tl  ds+  c>  (53) 

and  by  evaluating  the  integral  for  real  values  of  $ and  using  equation 
(1$),  which  is  valid  for  this  case 


t:=  Ca[^c  X*  fyff]*  C) 


(5U) 


* 


17. 


For  S c + / 


and  for  j*_  _/ 


y-  = — ^vi  607  ¥•  tt (. 


therefore 


^ ~ ~T r IF  ) + £/ — c^7  -a 
c0  = - Fr L 


(55) 


and 


Thus  the  result 


C , - JZ„  us,  + (jr  -r  fyz)u 


is  valid  for  1*1  < / * Equation  57  may  be  rewritten  as 

~ Xv\  ur  ■/  <-  h (S) 

where 

Aft)*  ir+fe-P/ir  «.«  i/>  J-  £«-«  ^ Inv  Kir* 

C.  Determination  of  the  Complex  Potential  Function: 

The  complex  potential  function  of  the  flow  in  the  upper  3 half 
plane  has  the  form 


(56) 

(57) 

(58) 

(59) 


X-  ft $ +"j&rr  -P'* 


(60) 


where  Q is  the  strength  of  the  total  sink  located  at  S~!  and  its  sign 
is  negative.  A source  of  equal  strength  is  located  at  infinity  giving 
X.  the  required  singularities.  Now 


-jij 


(61) 


becomes  zero  for  two  values  of  J . If  A and  B are  determined  so  that 
these  values  coincide  with  and  the  00  ential  X represents 


the  desired  flow.  Then 


zd/?c  + ~B  + -pry  — O 

Zfls  ■+  "B  t & t=T~  ~ 0 

A = (TTzi-Tr 

3=  ci-fyl-o 


(62) 


and 
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With  these  values,  equation  61  becomes 

- -2_  . (C  r , (Ah) 

' Z7T  (i-iXV-iJS-l)  {bh) 

Before  determining  the  shape  of  t..e  slot  a generalization  of  this  result 
is  in  order. 


J*3L 
d S 


The  preceding  investigation  has  been  based  upon  a consideration  of 
the  flow  near  the  trailing  edge  only.  Without  a suction  slot,  the  flow 
is  identical  to  the  stagnation  point  flow  represented  by  Figure  20.  This 
is  an  approximation  for  the  flow  at  the  trailing  edge  of  a profile  and 
no  influence  due  to  change  in  angle  of  attack  is  provided  for.  For  larger 
angles  of  attack  where  this  influence  has  to  be  taken  into  account,  a stag- 
nation point  can  be  added  to  the  flow  on  the  lower  surface  so  that  the 
flow  in  Figure  -20  would  become  that  of  figure  -21.  There  must  be  added 
a third  order  term  of  S to  the  right  side  of  the  complex  potential  func- 
tion X-  °?  equation  60.  If  So  is  the  S value  corresponding  to  » 
one  obtains  instead  of  equation  6ii  the  generalization 


Of  course  the  hodograph  is  more  complicated  in  this  case,  however, 
the  considerations  of  the  preceding  section  can  be  generalized  for  this 
also.  As  long  as  the  angle  of  attack  is  not  extremely  high  and  the 
forward  stagnation  point  of  the  profile  is  s till  situated  in  the  vicinity 
of  the  leading  edge,  this  generalized  solution  will  not  be  necessary. 


D.  Determination  of  the  Shape  of  the  Slot* 
The  hodograph  relation 

yields 
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Combining  equations  6 It  and  66  results  in 

"df  r $r  e 

The  function  t ( S ) is  known  from  part  B of  this  section,  so  z can 
be  found  as  a function  of  X from  the  above  equation. 

Substituting  into  equation  67  the  value  of  t as  obtained  by  equatidn 
$2  yields 

d-fk  _ O C-  JT  rs  - t + ✓**■/  If*-/  / 
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for  j5/~  I . In  the  vicinity  of  _J*=  / this  function  behaves  according  to 

Ql. 


ZTT  ur;  X - / 

Passing  over  S 3 / within  the  upper  J half -plane  from  a real  value 
*5  y / to  a real  value  / the  function  z is  augmented  by 


£ TTur, 

The  slot  width  is  then 
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£ = (69) 

The  same  expression  for  6"  results  if  one  uses  equation  57  for  t and  passes 
over  r / from  a value  3"  < / to  a value  ~$  > / , 

For  the  computation  of  the  shape  of  the  slot  and  the  velocity  distri- 
bution along  the  slot  walls,  the  following  expressions  can  be  used  which 
result  from  equations  67  with  5?  or  57s 

if  151=/ 
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and  if  / J/  = / 

X^^Ys^ih)  k°S  ',(V-CS'"  hM) 

where 
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The  function  z - -fft)  for  the  above  equations  has  to  be  determined  by 
graphical  or  numerical  integration,  the  z-values  of  the  curved  wall  of 
the  slot  corresponds  to  real  values  of  with  / J /£  / . 

The  velocity  distribution  along  the  walls  follows  from  equation  66 
which  yields  together  with  equations  52  and  57  for  Isis/ 


and  for  / J/~  / 


u^T 


(73) 


Figure  22  shows  a suction  slot  commuted  lor  the  trailing  edge  of  the 
NACA  23015  profile.  The  slot  width  £ has  been  chosen  equal  to  unity. 
Thus  */  from  equation  69.  The  angle  0 — 17.2°  and 

the  parameters  C and  S have  been  chosen^ 

c = 5 = -2  . 

Their  choice  determines  the  length  of  the  lower  lip  of  the  slot  and  the 
location  of  the  stagnation  point  <5*  . For  the  case  C~  S the  stagnation 
point  has  its  rearmost  position  and  is  ready  to  leave  the  profile  with 
an  increased  amount  of  suction. 


III.  DESCRIPTION  CF  EXPERIMENTAL  APPARATUS 


A.  Wind  Tunnel  and  Allied  Apparatus t 

All  of  the  quantative  measurements  of  lift,  drag,  pitching  moment 
and  pressure  readings  were  taken  in  the  Princeton  University  3^'  * 5'  sub- 
sonic wind  tunnel.  In  its  original  configuration  this  tunnel  was  three 
dimensional,  however  in  this  program  the  tunnel  was  converted  to  the  two 
dimensional  by  the  construction  of  a new  contraction  section,  test  section 
and  diffuser.  It  was  the  construction  of  the  diffuser  that  presented  the 
first  real  problem  in  that  flew  separation  from  the  walls  could  not  be 
avoided  and  the  location  of  the  initial  separation  point  could  not  be 
maintained  constant.  This  resulted  in  a tunnel  with  a quite  low  energy  ratio 
(which  could  be  tolerated  in  this  investigation)  and  with  severe  velocity 
fluxuations  in  the  test  section.  Of  course  the  test  velocity  varation  was 
unacceptable.  The  first  attempt  at  a "cure"  was  the  splitter  configuration 
as  shown  in  figure  23.  This  showed  a very  definite  improvement  in  the 
steadiness  of  the  flow  in  the  test  section.  However  the  problem  was  still 
largely  unsolved.  Further  attempts  with  slats  to  direct  the  diffuser  flow  to 
the  walls  and  spoilers  to  fix  the  position  of  the  stall  proved  unsatisfactory 

In  order  to  bring  the  operating  characteristics  of  the  two  dimensional 
test  section  up  to  an  acceptable  level  in  as  little  time  as  possible  two 
approaches  were  made  simultaneously.  They  were:  (1)  the  design  of  vortex 
generators  fran  data  supplied  by  the  United  Aircraft  Corp.  and  (2)  the 
design,  construction  and  installation  of  a new  low  angle  diffuser  equipped 
with  cusps  at  the  down  stream  end.  The  vortex  generators  were  finished 
last  but  not  installed  by  the  time  the  new  diffuser  with  cusps  was  completed. 
This  arrangement  was  found  to  be  wholly  satisfactory  and  the  tunnel  was 
operated  in  this  configuration  throughout  the  test  program.  (See  Figure  2h) 

The  wind  tunnel  balance  is  a six  component, oil  dampened,  pneumatic 
system  indicating  forces  and  moments  directly  on  large  differential  pressure 
gages.  For  this  investigation,  however,  only  the  components  of  lift,  drag 
and  pitching  moment  had  significance.  Drag  is  discussed  in  a later  section. 


T.  Reference  3 


however,  it  was  found  desirable  to  use  a wake  rake  in  order  to  determine 
drags  by  momentum  considerations.  This  rake  is  of  the  differentiating  type 
and  is  located  in  the  test  section  approximately  one  and  one  half  chord 
lengths  behind  the  trailing  edge  of  the  model.  The  total  head  and  static 
tubes  in  the  rake  were  connected  to  a multiple  manometer  equipped  with 
lights  and  a graflex  camera.  Thus  all  wake  drag  data  and  all  other  pressure 
data  were  taken  by  means  of  photography. 

B.  Smoke  Tunnel; 

Early  in  the  investigation  of  trailing  edge  suction  it  was 
realized  that  a number  of  unknown  qualities  of  the  flow,  associated  with 
this  type  of  high  lift  device,  might  be  revealed  by  a smoke  study.  This 
realization  eventually  led  to  the  purchase  of  a 3'x5'  smoke  tunnel  designed 
by  Dr.  A.  Lippisch  and  built  by  Collins  Radio  Corp.  of  Cedar  Rapids,  Iowa. 

It  was  felt  that  there  would  be  a considerable  time  delay  in  the  delivery 
of  such  an  apparatus  and  the  decision  was  made  to  build  a small  smoke  tunnel"1- 
in  order  to  solve  some  of  the  more  immediate  problems  of  circulation  control 
systems. 

The  smoke  tunnel  as  used  in  this  study  is  purely  a means  for  obtaining 
qualitative  information  since  the  flow  velocities  were  so  small  that  neither 
force  or  pressure  reading  could  be  taken  to  indicate  lift  or  drag.  Except 
for  this  the  tunnel  is  well  instrumented  with  a means  for  determining  flow 
velocities  and  suction  quantities.  This  instrumentation  offered  a basis  for 
accurate  comparison  of  the  flows  studied. 

C.  Models; 

The  primary  model  used  in  this  study  of  trailing  edge  suction  is 
a NACA  23015  profile  of  rectangular  planform  complete  with  end  plates. 

This  model  is  hollow  and  is  constructed  entirely  of  brass.  It  is  equipped 
with  22  chordwise  pressure  taps  and  with  nozzle  like  ends  outboard  of  the 
end  plates  through  which  the  air  removed  via  the  suction  slot  is  withdrawn. 

This  model  has  an  18"  span  between  the  end  plates  and  a base  chord  of  8", 
however  the  after  20%  of  the  chord  was  made  removable  for  the  attachment  of 
trailing  edge  suction  slots  (see  figure  25  & 26).  A total  of  ten  different 
combinations  of  upper  and  lower  suction  slots  were  attached  to  this  model 
and  completely  tested.  In  addition  a series  of  experiments  were  made  with  a 
variable  trailing  edge  also  made  of  brass  and  adjustable  between  wind  tunnel 
runs.  In  addition  to  the  brass  23015  airfoil  two  16  inch  chord  wood  models 
(see  figure  27)  were  built  and  tested  both  by  force  and  pressure  measurements. 

A 15  percent  thick  elliptical  model  was  also  built  to  obtain  preliminary 
information  on  the  behavior  of  the  bent  up  lower  lip  as  finally  used  ijn 
configuration  l*-f.  The  data  taken  with  this  model  is  reported  in  Reference  7. 
Other  models  used  in  the  study  were  5 inch  chord  2 inch  span  profiles  built 
for  testing  in  the  smoke  tunnel.  These  were  geometrically  similar  in 
trailing  edge  configuration  to  the  pressure  distribution  models  (see  figure  28). 
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D.  Suction  Systems 

For  the  wind  tunnel  tests  of  all  models,  except  those  built  for 
the  smoke  tunnel,  suction  was  obtained  by  using  an  auxilliary  stage 
compressor  taken  from  an  Allison  V-1710  engine.  This  compressor,  driven  by 
a UO  HF  D.C.  motor  was  located  outside  the  wind  tunnel  and  air  was  ducted 
from  the  model  b>  the  system  shown  in  figure  1$,  In  order  to  avoid  affecting 
the  balance  readings  by  introducing  loads  from  the  air  ducting^ a special 
coupling  consisting  of  two  ball  and  socket  joints,  connected  by  piping 
containing  a slip  joint^was  employed.  Flexible  piping  was  used  to  connect 
this  coupling  to  the  ducts.  Properly  located  (according  to  the  ASME  power 
test  code)  in  the  duct  is  a box  which  permits  insertion  of  several  different 
size  orifice  plates  for  measuring  the  flow  or  amount  of  suction.  Control  of 
the  amount  of  suction  is  obtained  by  varying  the  speed  of  the  D.  C.  motor 
and/or  by  changin£  the  vaJve  position  of  the  hydraulic  coupling  between  the 
power  input  to  the  compressor  and  the  compressor  wheel. 


I/.  METHODS  CF  DATA  HEDICTICN 


This  section  will  deal  with  the  methods  and  techniques  used  to  reduce 
the  raw  wind  tunnel  data  to  coefficient  form,  and  the  corrections  applied  to 
this  data  in  order  to  make  a clear  and  accurate  presentation. 

A.  Force  Measurements 


All  force  measurements  (lift,  drag,  and  pitching  moment)  were 
reduced  to  coefficient  form  on  the  basis  of  the  actual  chord  of  the  model 
tested  and  not  upon  the  base  chord,  ^Similarly  the  pitching  moment  coefficient 
is  based  upon  the  quarter  chord  length  of  the  actual  chord  and  not  the  base 
chord  of  the  profile. 

Wind  tunnel  test  section  velocities  were  calibrated  with  a pitot-static 
tube  for  each  model  tested.  This  was  accomplished  by  setting  the  model  at 
zero  degrees  angle  of  attack  with  the  pitot  tube,  located  in  the  middle  of 
the  tunnel  over  the  model,  mid  way  between  the  model  and  the  top  of  the  test 
section.  A calibration  was  made  between  the  pitot-static  readings  and 
the  difference  in  static  pressure  between  the  settling  chamber  and  contraction 
section.  These  latter  pressures  were  obtained  by  means  of  wall  taps  in  the 
respective  sections.  This  yielded  a calibration  curve  for  g.  plotted  against 
tunnel  differential  pressure  and  all  subsequent  speeds  were  measured  by  this 
AT  reading.  The  slope  of  this  calibration  c urve  together  with  data  taken 
during  the  tests  yielded  a simple  method  for  reduction  of  force  measurements 
to  coefficient  form.  The  wind  tunnel  balance  was  initially  calibrated^-  and 
the  results  plotted  so  that  each  of  the  forces  (i.e.  lift,  drag,  side  force) 
and  the  moments  (pitch,  ro-1,  yaw)  appear  as  a straight  line  with  force  or 
moment  8.  . ^age  reading.  Thus  the  slope  of  each  of  these  balance  calibration 
curves  multiplied  by  the  difference  in  wind  off  ano  wind  on  gage  readings 
yields  the  force  or  moment  directly.  Therefore  the  expression  for  lift 
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coefficient 


can  be  changed  to 
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■ the  difference  between  wind  on  and  wind  off 
lift  gage  readings 

* slope  of  lift  calibration  curve 

• wind  tunnel  differential  pressure 


/fe  * slope  of  j calibration  curve 
% 

Kl.  - ■ overall  lift  constant 

In  a similar  manner  the  expression  for  drag  and  pitching  moment  coefficients 
are  derived^resulting  in 
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It  is,  of  course,  of  great  importance  to  properly  evaluate  the  tares 
for  any  wind  tunnel  test.  The  force  model  was,  in  this  installation, 
mounted  on  end  plates  having  a diameter  somewhat  greater  than  the  chord. 
Balance  readings  indicated  excessively  high  drag  coefficients  but  quite 
normal  lift  coefficients^-  in  good  agreement  with  the  published  NACA  results. 
These  lift  values  were  obtained  with  the  conventional  sharp  trailing  edge 
with,  of  course,  no  suction.  Other  tests  performed  in  the  same  wind  tunnel, 
which  were  part  of  a project  being  carried  out  concurrently  with  the  subject 
of  this  report,  indicated  that  leakage  into  the  wind  tunnel  at  the  model 
mounts  had  a detrimental  effect.  This  effect  was  not  noticed  with  the 
circulation  control  force  model  at  zero  suction  coefficient.  However,  it 
seemed  reasonable  to  expect  a greater  leakage  where  suction  w«s  applied  .. 
with  the  possibility  of  introducing  some  error.  Means  were  not  available 
at  the  time  to  evaluate  properly  the  effect  of  this  leakage.  However, 
partial  pressure  distributions  taken  about  the  model  with  tue  end  plates 

II  This  is  explained  in  more  detail  in  section'  \/-A  of  the  discussion. 


aed  sx’»  *ee*ec  ladies  tax  ti*  lift  ie#*its  tn  ecELServative  by  a fa* 
percent.  Ik*  ccntlnaiaf  sacperiamsits  wile.  this  xcewi,  mi  eh  centimes  flap* 
»it£  trailing  edge  macticc,  will  inclade  a stody  to  evaluate  Use  teres  with 
suction.  As  previously  sectioned  the  balance  dreg  readings  ware  exetsslnlf 
cign.  As  so  expedient  solstice  the  drags  were  deiemiDec  by  sears  of  pressure 
surveying  uc  wake  . 

Slue*  tha  force  aets.reater.ts  were  takes  with  the  aodel  moon  ted  at  the 
SOt  base  coord  point  the  pitching  moment  coefficients  were,  cf  course,  about 
this  point.  These  data  were,  however,  ecu  verted  to  the  2$t  actual  chord 
point  Cf  each  configuration  by  scans  of  the  following  expression,  the 
signif icance  of  wnich  is  easily  understood  by  reference  tc  the  accompanlng 
diagram. 
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The  tars, base  chord,  refers  to  the  original  chord  of  the  unmodified  NACA  23015 
profile  while  the  actual  chord  refers  to  the  profile  with  the  sharp  trailing 
edge  removed  and  suction  slots  installed. 


B.  Preesure  Measuresente 


1.  Pressure  Distributions: 


The  experience  gained  in  testing  and  analysing  the  results  of 
the  force  measurement  model  answered  seme  questions  about  trailing  adge  suction 
but  created  many  new  ones.  The  question  of  the  type  of  pressure  distribution 
in  the  vicinity  of  the  slot  and  the  effect  of  the  increased  circulation 
( preseurewisa)  led  to  the  construction  and  testing  of  two  16”  chord  models 
equipped  with  pressure  taps.  The  results  of  these  tests  are  discussed  in 
section  VI. 


These  pressure  distributions  were  plotted  in  a direct  manner  (Fig.  30 
shows  a typical  plot  for  trailing  edge  1-b) . Since  the  data  ware  taken  by 
photographing  a multiple  mancme ter,  the  manometer  readings  themselves  were 
plotted  versus  the  chord  length  (In  inches).  This  yielded  an  area 
proportional  to  lift.  Mathanical  integration  of  the  area  ancloead  by  the 
curve  was  accomplished  by  use  of  a planimeter.  This  area  multiplied  by  a 
constant  relating  square  inches  under  the  curve  to  pounds  of  lift  resulted 
in  the  determination  of  tha  lift  coefficient.  When  the  preesure  distribution 
ia  plotted  in  this  manner  one  must  be  careful  to  establish  the  ordinate  of 
the  leading  adge  stagnation  point.  This  is  accomplished  by  plotting  a 
straight  horisontal  line  having  tha  ordinate  of  tha  teat  section  static 
preesure.  The  algebraic  addition  of  tunnel  | to  this  static  pressure 


1.  tt»ia  le  discussed  more  fully  in  tha  next  section. 


yields  the  ordinate  of  the  stagnation  point*  For  the  example  shown  in 
Figure  30  the  tunnel  static  pressure  was  -6*70  in.  HgO  (gage)  and  tunnel 
9 was  It. 31  in.  H2O.  The  stagnation  point,  therefore,  is  represented  by 
the  ordinate  2*39  in.  of  HpO.  Thus  the  stagnation  point  corresponds  to  a 
>Cp  of  *1.0  and  the  horizontal  line  representing  test  section  static  pressure 
corresponds  to  >Cp  ■ 0.  This  method,  rather  than  the  more  conventional  ^ 
transfer  of  the  data  to  pressure  coefficients,  was  used  because  it  eliminated 
a great  deal  of  routine  computation.  Since  the  pressure  distributions  made 
numbered  in  the  hundreds  a significant  amount  of  time  was  saved. 

2.  Equivalent  Power  Drag  Coefficient 

In  order  to  compare  the  slot  efficiencies  of  the  several 
configurations  the  following  procedure  was  undertaken*  A kiel  tube 
located  within  the  model  was  used  to  determine  the  total  pressure  on  the 
down  stream  side  of  the  slot.  The  external  total  pressure  was  determined  by 
a pitot  tube  placed  in  the  free  stream.  Thus  a pressure  coefficient, 
was  evaluated,  for  a range  of  <s,  for  each  trailing  edge  slot.  The  drop 
in  total  pressure  across  the  slot  is  a function  of  its  drag  and  hence  the 
power  requirement  for  a given  ^ C a • An  equivalent  power  drag  coefficient 
Xjg  was  evaluated  by  means  of  the  following  expression^ 

= OCf~0  • 

3.  Profile  Drag  Coefficient  by  Wake  Survey 

The  profile  drag  coefficients  presented  in  this  report  were 
obtained  by  a pressure  survey  of  the  wake 2.  A differentiating  wake  rake 
located  in  the  wake  of  the  model  was  used  in  combination  with  a multiple 
manometer.  Manometer  readings  were  taken  by  means  of  photography  yielding 
an  instantaneous  picture  of  the  reduced  total  head  within  the  wake.  The 
method  of  reduction  of  these  data  follows  the  procedure  outlined  in  Reference 
6 which  results  in  the  expression  for  drag  coefficient. 
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■ wake  width 

■ profile  chord 

• free  stream  dynamic  pressure 

* dynamic  pressure  within  the  wake 

■ incremental  distance  between  the 
wake  dynamic  pressure  readings. 
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The  shape  of  the  wake  pressures  when  seen  on  a multiple  manometer  is 
best  described  as  a bell  shape  curve  which  tends  to  spread  out  and  flatten 
with  increasing  angle  of  attack.  This  occurs  in  a manner  so  that  the  area 
under  this  curve  constantly  increases  but  at  a slower  rate  than  the  wake 
width  increases. 

With  increasing  suction  the  wake  width  and  the  area  under  the  bell 
curve  both  decrease  but  the  decrease  in  wake  width  is  the  more  pronounced. 

Early  in  the  circulation  control  study  the  procedure  for  solving  the 
wake  drag  equation  utilized  the  mechanical  integration  of  the  integral 
Sidy  by  means  of  a planimeter.  The  results  of  this  method  were  somewhat 
disappointing  in  that  a great  deal  of  scatter  of  the  test  points  was 
encountered.  In  order  to  arrive  at  a better  way  of  evaluating  the  integral* 
a plot  of  x * y^  was  carefully  made  between  limits  which  yielded  an  area  of 
approximately  the  same  magnitude  as  the  experimental  area.  A mathematical 
integration  of  this  area  was  made  and  used  as  the  standard  in  a comparison 
with  several  other  methods  of  integration.  These  other  methods  included: 

1.  Planimeter,  percentage  error  1 .2% 

2.  Trapezoidal  rule,  percentage  error  0.72% 

3.  Simpson's  rule,  percentage  error  Q.V?% 

The  percentage  error  is  based  upon  a zero  percentage  error  for  the 
mathematical  integration.  While  none  of  the  errors  look  excessive  at  first 
glance  nevertheless  an  inspection  of  the  equation  7 h shows  that  the  drag 
coefficient  is  a small  difference  of  two  relatively  large  numbers.  Thus 
any  error  involved  in  evaluating  the  integral  is  magnified  many  times  in 
the  drag  coefficient.  Accordingly,  a procedure  was  set  up  whereby  / 3 d1 
could  be  evaluated  automatically  on  an  I.B.M.  computing  machine  set  to  use 
Simpson's  Sule.  In  fact  it  was  found  to  be  economically  sound  to  feed  all 
of  the  wake  information  into  the  machine  and  receive  the  wake  drag 
coefficients  directly. 

C.  Corrections  hade  to  the  Wake  Drag  Coefficient: 

One  of  the  interesting  characteristics  of  trailing  edge  suction  is 
that  it  tends  to  remove  the  effect  of  friction.  In  other  words  the  real  flow 
tends  to  approach  potential  flow  in  some  respects.  Thus  as  the  suction 
quantity  is  increased, more  and  more  of  the  turbulent  wake  is  swallowed  into 
the  suction  slot,  in  fact  under  some  conditions  of  high  jC*  it  almost 
entirely  disappears  and  the  wake  drag  coefficient  approaches  zero.  This  can 
only  be  the  profile  drag  coefficient  under  one  condition  and  that  is  the 
removed  air  must  be  returned  to  the  free  stream  so  that  there  is  no  change 
in  momentum  between  the  slot  and  the  component  in  the  downstream  direction 
of  the  discharged  air. 

If  this  is  not  the  condition  then  there  must  be  a correction  made  to 
the  wake  drag  to  yield  the  profile  drag  coefficient.  The  development  of  this 
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correction  is  as  follows:  Figure  31  shows  a schematic  sketch  of  an  airfoil 

equipped  with  a trailing  edge  suction  slot.  Streamline  ^ having  a velocity 
u/i  at  station  1 stagnates  at  station  2 from  2 it  accelerates  thru  the  slot 
then  turns  at  3 and  flows  in  a spsnwise  direction.  At  point  3 the  velocity 
component  in  the  direction  of  the  free  stream  is  zero. 

Therefore: 


At  station  1 the  velocity  ■ u/© 

At  station  2 the  velocity  ■ 0 

At  station  3 the  velocity  ■ 0 

The  momentum  change  between  stations  1 and 

The  momentum  change  between  stations  2 and  3*0 

Therefore  the  increased  drag  due  to  momentum  change  of  the  removed 
air  is 


-/*Q 

and  the  increased  drag  coefficient  is 

sCd  = rAs  = * Z/C« 


D.  Measurement  and  Seduction  of  Suction  Quantity  Data: 

The  procedure  by  which  flow  quantities  were  measured  and  the  data 
reduced  is  that  outlined  in  the  A.S.M.E.  Power  Test  Code;  The  sharp  edged 
orifice  was  selected  as  the  primary  element  because  of  its  reliability  and 
low  cost.  It  is  true  that  such  a metering  device  is  less  efficient  than  a 
nozzle  but  an  excess  of  power  was  available  and  it  appeared  more  economical 
to  spend  power  than  to  suffer  the  cost  of  a nozzle. 

This  method  of  determining  the  flow  quantity  is  in  essence  based  upon 
the  fundamental  flow  equation 


(75) 

where  £*P  is  the  pressure  drop  across  the  sharp  edge  orifice.  For  the 
test  installation  vena  - contracta  pressure  taps  were  used  to  measure  this 
pressure  drop.  The  actual  flow  equation  used  by  reference  7 relates  mass 
flow,  uj  i to  the  parameters  of  equation  75*  In  order  to  obtain  Q 
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(ft?/sec.,  tunnel  density  was  obtained  by  measuring  settling  chamber  static 
pressure  and  temperature  and  dividing  this  result  into  the  mass  flow,  c o . 

This  yielded  the  suction  quantity  based  upon  the  tunnel  conditions  rather 
than  upon  those  of  the  duct.  The  A.S.M.E.  procedure  was  carefully  followed 
for  a number  of  tests  including  various  angles  of  attack  and  several  different 
slot  shapes  and  dimensions.  It  was  then  found  that  a calibration  curve 
could  be  constructed  relating  *47  across  the  orifice  directly  to  mass  flow 
thru  the  orifice  and  that  the  error  involved  in  the  curve  was  well  within 
the  experimental  error  of  obtaining  the  data.  Such  a calibration  curve 
indicates  that  the  density  variation  with  A~P  , at  the  orifice,  was 
approximately  the  same  regardless  of  trailing  edge  geometry.  With  this 
calibration  established  it  was  then  possible  to  greatly  speed  up  the  testing 
program  and  the  data  reduction.  The  suction  quantity  was  non-dimens ionalized 
into  coefficient  form  by  means  of  the  relationship, 

^ &_ 

S^q  — cji  s 

Here,  as  in  the  case  of  all  other  coefficients, & q was  computed  on 
the  basis  of  the  actual  area  tested  and  not  upon  the  base  area  of  the  profile. 

V.  DESCRIPTION  CF  EXPERIMENTS 

A.  Calibration  of  the  Wind  Tunnel 

As  previously  mentioned  the  Princeton  3£,x5'  subsonic  wind  tunnel 
was  modified  for  this  experimental  program  by  the  installation  of  an 
auxilliary,  two  dimensional  18"  throat.  After  the  initial  velocity  variation 
difficulties  had  been  overcome  by  the  installation  of  the  cusped  diffuser  it 
became  necessary  to  examine  the  character  of  the  flow  in  the  new  test  section 
in  order  to  make  proper  corrections  to  the  test  data.  This  investigation 
was  made  in  four  parts: 

1.  Velocity  survey  of  the  test  section. 

2.  law  head  survey  of  the  region  of  the  test  section  where  the  model 

would  be  installed. 

3.  Determination  of  the  turbulence  of  the  wind  tunnel. 

ii.  Measurement  of  the  boundary  layer  on  the  side  walls  of  the  test 
section  in  the  vicinity  of  the  model. 

The  velocity  survey  of  the  test  section  was  made  by  probing  horizontally 

with  a pitot  static  tube  at  different  elevations  f rcm  the  bottom  of  the  t est 

section  to  the  top.  All  of  the  points  tested  were  outside  the  boundary  layer 
end  the  results  were  considered  good  since  the  maximum  and  minimum  values 
of  3 were  within  1.0J(  of  the  average  ^ • 

In  order  to  determine  the  inclination  of  the  flow  in  the  test  section 
a yaw  head  survey  was  made  across  the  mid  point  of  the  test  section.  The 
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average  deviation  of  the  flow  from  the  horizontal  was  approximately  one 
half  a degree  and  the  inclination  was  upward.  However,  with  the  model 
installed  at  -1°  to  the  horizontal  the  inclination  of  the  flow  was 
approximately  zero  degrees.  Since  the  angle  of  zero  lift  for  the  base  23015 
profile  is  -1°  and  repeated  tests  with  this  model  gave  excellent  agreement 
with  the  NACA  results  it  was  decided  that  the  angularity  correction  was 
essentially  zero. 

The  turbulence  factor^-  was  determined  by  use  of  the  conventional 
turbulence  sphere  fitted  with  fcur  static  pressure  taps  on  the  downstream 
side  22°  frc*u  the  horizontal.  The  results  indicated  a critical  Reynolds 
number  of  357.200  which  yielded  a turbulence  factor  of  1.078.  The  reason 
for  this  rather  low  turbulence  is  probably  the  large  contraction  ratio 
(seventeen)  of  the  auxiliary  throat  and  the  effect  of  the  cusps  in  permitting 
a small  angle  diffuser. 

Information  about  the  thickness  and  character  of  the  side  wall  boundary 
layer  was  necessary  to  determine  what  corrections,  if  any,  would  be 
necessary  to  obtain  the  mean  dynamic  pressure  in  the  test  section.  This 
was  accomplished  by  means  of  a total  head  survey  of  the  boundary  layer  at 
the  location  of  the  model  in  the  test  section.  The  assumption  of  constant 
static  pressure  across  the  boundary  layer  permits  the  relating  of  total 
head  within  the  boundary  layer  to  the  dynamic  pressure.  The  investigation 
revealed  a boundary  layer  at  the  model  location  of  approximately  one  inch 
in  thickness.  However,  this  was  a turbulent  condition  and  the  difference 
between  the  calibrated  free  stream  and  the  average  test  section  based 
upon  this  turbulent  layer  was  of  an  order  of  magnitude  that  could  be 
neglected.  Once  again  the  justification  for  a lack  of  correction  is  the 
excellent  agreement  between  results  with  this  tunnel  and  those  published 
by  the  NACA  not  only  for  the  base  23015  profile  but  for  other  profiles  more 
recently  tested. 

B.  Calibration  of  the  Blower; 

As  mentioned  earlier  in  this  report  the  duct  from  the  wind  tunnel 
to  the  compressor  contained  an  apparatus  permitting  the  rapid  change  of 
orifice  plates.  Three  different  orifices  were  originally  used  to  measure 
the  flow  over  a wide  range  of  flow  quantities.  These  were  1",  2"  and  1*.9W 
in  diameter  with  & ratios  of  0.11*3>  0.286  and  0.70  respectively.  After 
testing  several  t railing  edge  configurations  thru  wide  range  of  angles  of 
attack  it  was  possible  to  construct  the  calibration  curve  discussed  in  the 
section  on  data  reduction.  As  a check  on  the  validity  of  the  calibration 
a duct  was  attached  to  the  exhaust  part  of  the  blower  and  the  flow  in  this 
duct  was  measured  with  an  anemometer.  Pressures  and  temperatures  were 
taken  simultaneously  with  tne  velocity  readings.  Prom  these  data  mass  flow 
was  computed  which  gave  excellent  agreement  with  the  calibration  for  both 
the  2”  and  1*. 9"  orifices.  Data  taken,  however,  for  the  1-inch  orifice  was 
not  in  good  agreement  with  its  calibration.  Several  explanations  for  this 
condition  are  possible  but  an  investigation  would  have  been  necessary  to 
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substantiate  any  one  of  them.  As  an  expedient  solution  it  was  decided  to 
discard  the  1-inch  orifice  altogether  especially  since  the  flow  coefficients 
were  not  as  reliable  for  this  /Q  ratio  as  for  the  larger  values. 

C.  Determination  of  Spanwise  Variation  of  Slot  Velocity! 

One  of  the  first  tests  made  with  the  force  model  was  the 
determination  of  the  spanwise  velocity  distribution  in  the  suction  slot. 

This  was  accomplished  by  the  installation  of  static  pressure  taps  on  the 
lower  lip  of  the  slot.  It  was  found  that  the  static  pressures  across  the 
slot  were  essentially  the  same.  With  the  assumption  of  constant  total 
pressure  there  was,  then,  practically  no  spanwise  variation  of  slot  velocity. 
This  was,  no  doubt,  due  to  the  physical r elationship  of  slot  size  to  the 
large  hollow  interior  of  the  model  and  the  large  ducts.  In  other  words,  the 
greatest  pressure  arop  in  the  system  occurred  across  the  suction  slot  so  this 
in  effect  became  a severe  throttling  device^ thus  tending  to  automatically 
distribute  the  flow.  This  large  pressure  drop  in  same  later  tests,  with  very 
small  slot  thicknesses,  caused  severe  bending  and  vibration  of  the  lower  slot 
lip  which  was  solved  by  the  installation  of  spacers  in  the  slot. 

D.  Model  Testing: 

In  general  the  various  configurations  of  the  force  model  were 
tested  in  the  following  manner.  With  the  tunnel  doors  closed  and  with  zero 
velocity  balance  readings  were  taken  of  lift,  drag  and  pitching  moment  for 
each  of  the  several  angles  of  attack.  The  model  was  then  returned  to  the 
desired  angle  and  the  wind  tunnel  brought  up  to  test  velocity.  With  a 
suction  quantity  of  zero,reaaings  of  lift,  drag,  and  pitching  moment  were 
again  taken.  This  process  was  repeated  for  approximately  fifteen  different 
suction  quantities  for  each  angle  of  attack  tested.  Concurrently,  the  wake 
pressures  were  photographed  on  the  multiple  manometer  for  each  suction 
quantity  for  each  angle.  Additional  data  taken  included  the  test  section 
and  settling  chamber  static  pressures,  the  settling  chamber  and  duct 
temperatures  and  the  absolute  and  differential  pressures  across  the  orifice 
plate  in  the  duct. 

A separate  series  of  tests  were  run  on  each  configuration  in  order  to 
measure  the  total  pressure  drop  across  the  slot.  This  was  done  so  the 
pressure  leads  from  the  model  to  the  manometer  could  not  influence  the  force 
measurements.  These  total  pressures  were  also  photographed  for  approximately 
twelve  different  suction  quantities  for  each  angle  of  attack  tested. 

A detailed  breakdown  of  the  configurations  tested  and  the  reasons  follow. 
Figure  32a,  b,  and  c shows  the  geometry  of  the  various  configurations  with 
the  exception  of  a variable  trailing  edge  which  will  be  subsequently 
discussed  in  this  section  and  in  section  VI.  Trailing  edges  1-a,  1-b,  and 
1-c  were  the  first  three  shapes  tested.  These  represented  the  range  of 
lower  lip  overhang  over  which  it  was  felt  a detailed  investigation  was 
warranted.  The  upper  lip  was  simply  formed  with  a single  radius,  for 
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although  it  was  felt  the  slot  shape  would  strongly  influence  the  results 
the  theory  of  such  a shape  as  presented  in  this  report  was  not  at  that  tine 
completely  developed,  lhe  test  results  for  these  first  three  models  indicated 
that  the  geometry  of  the  trailing  edge  strongly  influenced  the  maximum  lift 
of  the  section.  The  highest  lift  was  achieved  with  trailing  edge  1-b  which 
had  more  overhang  than  1-c  but  less  than  1-a  as  can  be  seen  in  Fig.  3?. 

Since  the  slot  thickness  varied  for  each  of  these  cases  it  was  thought  that 
there  might  be  an  optimum  ratio  of  overhang  to  slot  thickness.  Theory 
indicates  that  increases  as  the  slot  approaches  the  trailing  edge 

of  the  base  profile.  Consequently,  it  was  thought  that  in  the  case  of  1-c 
the  ratio  was  far  below  the  optimum.  As  a result  this  ratio  was  increased 
for  trailing  edge  2-b. 

Trailing  edge  2 is  the  "correct"  slot  shape  as  computed  by  the  method 
outlined  in  this  report.  This  was  first  tested  with  lower  lip  b and  the 
results  were  encouraging  in  as  much  as  the  at  low  angles  and  low 

suction  quantities  was  the  best  that  had  so  far  been  achieved.  The  maximum 
lifts,  however,  were  very  disappointing.  It  was  thought  that  flow  around 
the  lower  lip  might  account  for  this,  consequently  trailing  edges  d and  e 
were  designed  with  a longer  overhang  to  overcome  this  difficulty.  It  is 
interesting  that  again  the  maximum  lifts  achieved  with  the  section  were  with 
the  moderate  overhang  (trailing  edge  2-e). 

The  amount  of  air  that  could  be  sucked  through  the  slot  was,  in  this 
case,  small  compared  to  the  amount  for  the  No.  1 series  of  trailing  edges. 

It  was  thought  that  this  might  be  the  result  of  a flow  separation  at  the 
slot  inlet  so  trailing  edge  2 was  modified  internally  in  an  attempt  to 
effect  a cure.  This  resulted  in  trailing  edge  h. 

It  will  be  noted  in  Figure  32  that  trailing  edge  1*  has  an  internal 
cusp  to  aid  in  diffusing  the  flow  inside  the  model  and  to  prevent  any 
disturbance  from  working  upstream  around  the  inlet  contour.  The  bent  up 
lower  lip  of  configuration  ii-f  was  an  attempt  to  limit  the  stagnation  point 
travel  and  provide  easy  entrance  for  the  air  into  the  slot,  lhe  results 
were  mediocre  compared  to  1-b,  although  it  shows  a higher  maximum  lift  than 
2— e. 


Trailing  edge  5-e  is  another  application  of  the  cusp,  this  time  to 
provide  a more  efficient  inlet  contour  to  the  suction  slot.  The  results 
were  not  outstanding  either  in  maximum  lift  achieved  or  in  power  required 
to  move  the  air;  however,  recent  developments  indicate  that  the  cusp  in 
trailing  edge  5>  was  improperly  designed. 


The  two  configurations,  d-a  and  1-g,  were  tested  to  determine  the  effect 
of  extreme  shapes  such  as  the  sharp  upper  lip  in  trailing  edge  d-a.  The 
maximum  lift  for  this  profile  is  poor,  but  it  does  compare  more  favorably 
with  some  of  the  other  models  at  lower  's. 


The  experience  of  testing  these  models  revealed  that  there  can  be  a 
substantial  flow  around  the  lower  lip  into  the  slot  under  certain  conditions 
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of  angle  of  attack  and , this  forming  a free  stream  stagnation  point. 
Trailing  edge  1-g  was  designed  to  resist  this  tendency,  but  the  results 
. indicate  that  1-g  is  not  as  effective  as  trailing  edge  1-b  from  which  it  was 

formed. 

The  results  of  these  tests  and  the  general  theory  indicated  that  the 
location  of  the  slot  (relative  to  the  trailing  edge  of  the  base  airfcil)  and 
the  amount  of  overhang  had  substantial  effect  upon  the  achieved. 

In  order  to  experimentally  check  the  overhang  portion  of  this  hypothesis  a 
special  lower  lip  was  designed  and  constructed  which  was  adjustable  to  ten 
different  lengths.  This  lip  in  combination  with  upper  lip  No.  1 formed 
trailing  edge  1-v.  Figure  33  shows  the  geometry  for  all  ten  of  the 
adjustments.  Each  was  tested  for  lift,  drag,  and  pitching  moment  as  indicated 
by  balance  readings.  This  was  in  many  respects  the  most  interesting  experiment 
of  the  program  because  the  effect  of  systematically  changing  one  variable 
is  readily  understood. 

There  were  two  major  pressure  distribution  studies  made  during  this 
program.  These  involved  models  previously  described.  One  of  the  pressure 
models  was  geometrically  similar  to  configuration  1-b, the  other  was  a 
typical  trailing  edge  suction  model  designed  I-A.  In  each  case  measurements 
» were  made  following  the  procedure  described  for  the  force  model  since  one  of 

the  reasons  lor  these  tests  was  to  compare  the  results  of  lift  by  force 
measurement  and  pressure  distribution.  The  actual  testing  to  obtain  the 
local  static  pressures  around  the  profile  was  again  separate  from  the  force 
tests  so  that  pressure  leads  would  not  interfere  with  balance  readings. 

* All  the  pressures  were  indicated  upon  the  multiple  manometer  which  was 

photographed  for  several  suction  quantities  for  each  angle  of  attack  tested. 

A comparison  of  the  results  obtained  is  discussed  in  the  next  section. 


The  third  and  one  of  the  most  revealing  series  of  experiments  was  that 
of  studies  made  in  a smoke  tunnel^.  This  investigation  was  limited  to  the 
three  simplest  trailing  edges  tested  (1-a,  1-b,  and  1-c)  because  of  the  small 
scale  to  which  the  profiles  were  reduced  (*$  inch  chord).  The  smoke  tunnel 
results,  even  though  qualitative, are  difficult  to  compare  directly  with  the 
force  model  results  due  to  the  great  difference  in  Reynolds  Number.  It  was, 
however,  this  series  of  experiments  which  confirmed  the  free  stream 
stagnation  point  under  certain  conditions  of  and  angles  of  attack  and 
showed  that  this  condition  can  be  hastened  by  shortening  the  overhang 
(see  Figure  h) 


VI.  PRESENTATION  AND  ANALYSIS  CP  RESULTS 


% 


Detailed  plots  of  all  data  taken  have  been  prepared  although  the 
greater  part  of  these  are  presented  in  a separate  report1.  Individual 
plots  oi  sCjg  and>Cp  versus X $ for  each  angle 
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of  attack  have  been  made.  From  these  data  saui.sry  families  of  / Cx  versus 
angle  of  attack  at  several  values  of  have  been  plotted  for  each  trailing 
edge  configuration  and  these  are  presented  in  Figures  35  thru  1*6  representing 
the  lifts  as  obtained  by  force  measurements,  similar  families  are  shown  in 
Figures  1*7  and  lt6  but  in  this  case  the  lift  coefficients  were  obtained  by 
pressure  measurements.  A comparison  will  subsequently  be  drawn  between  these 
two  figures  and  Figures  1*5  and  1*6. 

Also  summarized  from  the  detailed  data  are  families  oZyCj  and 
versus  angle  of  attack  for  each  configuration.  These  a re  presenteu  in 
Figures  1*5  through  60. 

In  a similar  way  curves  representing  the  general  variation  of  jC  y>-  c/4 
with  angle  of  attack  ano  Xc^  are  presented  in  Figures  61  through  70. 

Figures  71  through  80  represent  the  variation  of^  (based  upon  total  head 
drop  across  the  slot)  with  angle  of  attack  for  sevfral  values  of  . Five 
additional  figures  complete  the  presentation  of  the  summarized  date.  These 
(Figures  81,  82,  83,  81* , and  85)  compare  the  characteristics  of  lift,  drag 
and  pitching  moment  of  the  various  trailing  edges  tested  and  figure  85 
shows  the  step  by  step  alteration  of  the  lift  coefficient  of  trailing  edge 
1-v  (variable)  as  the  overhang  is  systematically  increased. 

One  of  the  first  tests  made  with  the  force  model  was  with  a plain 
trailing  edge  in  order  to  compare  the  results  with  those  of  the  NACA  for 
the  same  profile.  Figure  31*  shows  this  comparison,  the  results  being 
identical  up  to  an  angle  of  attack  of  7°  where  the  great  difference  in 
Reynolds  Numbers  between  the  two  tests  begins  to  produce  its  effect.  It 
was  felt  that  this  test  was  an  excellent  overall  check  of  the  apparatus  and 
technique  used. 

Table  I compares  two  of  the  fundamental  characteristics  of  the  basic 
lift  curve  to  tnoae  of  the  lift  curve  for  the  several  trailing  edge 
configurations.  It  is  interesting  to  note  that  the  angle  of  zero  lift  for 
each  of  the  configurations  tested  with  the  force  model,  with  the  exception 
of  1-g,  is  identical  to  the  angle  of  zero  lift  for  the  base  profile. 

Trailing  edge  1-g  has  a severely  bent  down  lower  lip  which  would  lead  one 
to  believe  that  would  be  lower  than  for  the  other  cases)  however, 

the  flow  on  the  upper  surface  near  the  trailing  edge  is  probably  badly 
separated  due  to  both  the  bent  lower  lip  and  the  discontinuity  of  the 
upper  surface.  The  force  measurements  taken  with  the  pressure  distribution 
model  indicate  a slightly  lower  angle  of  zero  lift  than  for  the  force  model. 
This  may  be  due  to  small  differences  between  the  two  models,  perhaps  in 
the  curvature  of  the  mean  camber  line  especially  near  the  trailing  edge. 

The  angle  of  zero  lift  as  obtained  from  pressure  studies  is  slightly  higher 
than  that  obtained  from  force  measurements  with  the  same  model.  This  is 
due  to  the  fact  that  with  force  measurements  the  reaction  of  the  entire 
model  is  being  measured  but  a pressure  measurement  reveals  the  characteristics 
at  only  one  spanwise  station.  It  is  felt,  however,  that  the  agreement,  for 
this  caee,  is  quite  good.  An  examination  of  Table  I shows  that  the  slope 
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BASE  AIRFOIL  - NACA  23015 


TRAILING  EDGE 


FORCE 

hCQEL 


PRESS. 

iXDEL 


N.A.C.A. 

PLAIN 


A-a 

1-b 

1- c 

2- b 
2-d 
2-e 

l*-f 

5-e 

d-a 


vA-g 

ri-A 

I-B 
I-A 

I-B 


Force 

Force 

Pressure 

Pressure 


SLOPE  OF  LIFT  CURVE 

(d<f/do < ) 

0.108 

0.110 

0.113 

0.118 

0.121* 

O.lll* 

0.103 

0.123 

0.116 

0.120 

0.105 

0.109 

0.099 

0.106 

0.107 

0.106 


ANGLE  OF  ZERO  LIFT 
(DEG.) 

-1.0 

-1.0 

-1.0 

-1.0 

-1.0 

-1.0 

-1.0 

-1.0 

-1.0 

-1.0 

-1.0 

-0.5 

-1.3 

-1.5 

-1.2 

-1.0 


of  the  lift  curve  (^<0  varies  sar.ewhat  with  the  trailing  edge 
configuratioruhowever  there  does  not  seen,  to  be  enough  of  a trend  to  make 
deductions.  The  average  ( for  the  force  model  is  0.111*5  and 
for  the  pressure  model  0.102.  The  lormer  is  somewhat  higher  than  the  slope 
of  the  lift  curve  for  the  base  profile  as  given  by  the  NACA  while  the  latter 
is  slightly  lower.  The  end  plates,  which  exceeded  the  chord  of  the  force 
model,  did  not  completely  cover  the  ends  of  the  larger  pressure  distribution 
model.  It  is  possible  that  normal  leakage  around  the  end  plates  might 
simulate,  to  a degree,  three  dimensional  flow  which  might  account  for  the 
reduced  slope  of  the  lift  curve. 


Figures  35  thru  1*1*  representing  the  two  dimensional  lift  characteristics 
of  the  several  configurations  are  not  consistent  with  respect  to  maximum 

• This  is  due  to  sane  limiting  condition,  in  each  case,  of  the  amount 
of  suction  that  could  be  applied.  There  were  three  conditions  thus 
limiting  the  scope  of  the  tests.  The  first  was  encountered  with  trailing 
edges  1-c  and  2-b  and  this  was  the  establishment  of  a free  stream  stagnation 
point  after  which  the  lift  increments  went  in  the  negative  direction. 

Secondly,  severe  vibration  and  deflection  of  the  trailing  edge  was  encountered 
with  seme  of  the  configurations  but  this  was  overcane  by  the  use  of  spacers 
between  the  lips.  The  third  llmitationwas  the  size  of  the  slot,  the  greater 
the  slot  thickness  the  greater  th r ^Cq  • This  last  condition  limited  the 
majority  of  the  tests  and  it  represents  the  maximum  effort  of  the  system 
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for  each  trailing  edge  geometry.  Tfois  maximum  effort  should  not  be  confused 
with  a constant  power  condition)  merely  because  the  compressor  controls  were 
in  the  same  approximate  setting,  since  the  compressor  blades  operated  at 
different  efficiencies  for  each  slot  dimension.  In  view  of  this, one  would 
expect  that  the  highest  to  be  achieved  would  have  been  with  trailing 
edge  1-c  since  it  had  the  largest  slot  width.  This  case  was,  however, 
limited  by  a free  stream  stagnation  point  due  to  the  short  overhang  of  the 
lower  lip.  The  highest  ^9  was  reached  with  trailing  edge  1-b  which  had 
the  second  largest  slot  width  and  an  overhang  long  enough  to  prevent  the 
occurence  of  the  free  stream  stagnation  point.  This  discussion  has  so  far 
excluded  a comparison  with  the  g ecmetries  of  trailing  edge  1-v  which  is 
considered  to  be  a special  case.  Ibis  will  be  subsequently  examined. 

The  highest  lift  coefficient  measured  was  with  trailing  edge  1-b  and 
its  value  was  2.6 7.  This  was  only  partly  due  to  the  fact  that  1-b  achieved 
the  highest>£^  , since  examination  of  Figures  35,  36,  and  reveals  that 
this  trailing  edge  (1-b)  was  more  effective  st  comparable  >Cq  • s and  angles 
than  either  1-a  or  1-g  (it  should  be  remembered  that  1-g  is  1-b  with  the 
lower  lip  cambered).  This  indicates  that  is  influenced  by  the 

geometry  of  the  trailing  edge  suction  slot.  Figure  81  shows  a comparison 
of  the  different  configurations.  A study  of  this  figure  shows  that  the 
effectiveness  of  any  trailing  edge  depends  upon  the  angle  of  attack  and/Co 
For  example,  trailing  edge  5-e  with  the  external  cusp  has  the  highest  AA 
(0.83  at  +1©  angle  of  attack  for  a jCq  of  0.01  also  the  value  of  >£* 

(1.0)  for  this  suction  quantity  ana  angle  was  the  highest  tested.  For  the 
sam at  13°  angle  of  attack  trailing  edge  d-a  indicates  the  highest  lift 
coefficient  (1*5)  • It  is  seen  that  both  1—g  and  2— e exceed  the  above  value 
for  d-a  but  at  an  angle  of  11°.  Information  on  the  characteristics  of 
both  2-e  and  1-g  at  13°  angle  of  attack  is  not  availabla  However  the  lift 
summary  curves  of  Figures  35  thru  U*  indicate  that  these  higher  angles  for 
the  two  cases  mentioned  would  be  of  little  interest.  The  reason  that  the 
11  angle  of  attack  curves  of  trailing  edges  2-e  and  1-g  compare  so 
favorably  with  the  13°  curve  of  d-a  is  because  the  lift  coefficient  of 
d-a  at  13°  and  ■ .01  represents  the  partially  stalled  condition. 

Of  those  configurations  that  reached  a of  .02  the  maximum  lift  was 
achieved  by  trailing  edge  1-g  (Xjl  ■ 1.96)  at  11°.  This  is  again  due  to 
the  fact  that  1-g  had  not  stalled  at  11°  with  a XT*  - .02  while  1-b,  for 
example,  had  stalled  at  13°  with  the  same^C,  . Of  the  tests  during  which 
a Xtq  of  .01*  or  greater  was  obtained  trailing  edge  1-b  was  by  far  the  most 
effective.  Figure  35  indicates  that  the  reason  for  this  is  that  between 
* of  .02  and  .03  the  flow,  which  had  been  separated  at  13°  at  lower 
^C^'s,  becomes  attached.  Thus  the  effect  of  suction  is  not  only  an  upward 
displacement  of  the  lift  curve  but  is  also  an  extension  of  the  curve.  By 
such  behavior  trailing  edge  suction  at  the  higher  values  of  ^C9  simulates 
to  a degree  the  effect  of  both  slots  and  flaps. 

Figures  1*5  end  1*7  represent  the  lift  characteristics  of  trailing  edge 
I-A.  This  was,  as  previously  explained,  a 16"  chord  pressure  distribution 
model. However  force  tests  were  also  made  with  it  in  order  to  compare  the 
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results.  I-A  should  not  be  compared  to  1-a  since  the  trailing  edge 
proportions  are  substantially  different.  It  will  be  observed  in  figure 
u5  that  a distinct  change  in  slope  of  the  lift  curve  occurs  at  approximately 
5°  angle  of  attack.  This  is  shown  foiyCq  • o and  it  is  for  this  reason  force 
measurements  were  not  made  beyond  5 degrees;  however  pressure  distribution 
studies  were  made  up  to  13  degrees  angle  of  attack  and  the  shape  of  these 
curves  (Figure  Ii7)  are  more  normal.  It  is  seen  that  at  all  angles  and  c ' s 
the  results  from  the  pressure  investigation  are  substantially  greater  than 
those  of  the  force  measurements.  The  pressure  taps  were  located  at  the  mid- 
span  of  the  mcdel  and  were  therefore  less  influenced  by  the  separation 
which  began  at  the  side  walls  of  the  tunnel  and  worked  inward.  Of  course  as 
soon  as  this  separation  began  the  force  measurements  were  affected  and  this 
seems  to  have  started  at  about  5 degrees  angle  of  attack.  The  same 
explanation  holds  for  the  pressure  and  force  measurements  made  with 
configuration  I-B, the  results  of  which  are  shown  in  figures  1*6  and  1*8.  This 
trailing  edge  (I-B)  is  geometrically  similar  to  1-b  but  a comparison  of  the 
results  shows  1-b  to  be  far  superior  at  all  angles  of  attack  and  suction 
quantities.  Again  the  explanation  for  the  difference  between  the  8- inch 
chord  model  and  the  16 -inch  chord  seems  to  be  the  separation  which  begins 
at  the  wall  and  works  inward,this  condition  being  much  more  severe  for  the 
case  of  the  larger  chord. 

Figures  1*9  through  60  present  the  drag  characteristics  for  the  several 
trailing  edges.  Each  of  the  families  of  drag  curves  seem  to  have  the  same 
general  behavior  with  the  exception  of  those  of  1-c  and  2-b.  Figure  1*9  which 
represents  the  drags  of  trailing  edge  1-a  is  typical.  For /Cq  ■ o the 
vs.  oc  curve  is  rather  flat  until  an  angle  of  attack  of  5 degrees  has  been 
reacned.at  which  point  the  curve  abruptly  steepens.  As  the/?Q  is  increased 
the  whole  curve  is  displaced  downward  and  the  increase  of  drag  coefficient 
with  angle  of  attack  is  not  quite  so  steep.  It  can  be  seen  that  the  wake 
drag  coefficient  approaches  the  zero  drag  condition  at  the  higher 

of/C*  . As  explained  in  the  section  on  data  reduction  a value  of  2 jC<o 
was  added  to  the  wake  drag  for  momentum  considerations.  Curves  for  this 
are  also  plotted  with  the  wake  drag  coefficients.  It  should  be  realized 
that  the  * 2jC$  ) coefficients  are  valid  only  for  the  condition 

of  discharging  the  removed  free  stream  and  boundary  layer  air  in  a s panwise 
direction.  For  any  other  condition  the  drag  is  - £ drr\  where  jCdrr\ 

is  the  momentum  drag  coefficient.  This  may  be  increased  by  discharging 
the  air  forward  and  decreased  by  blowing  downstream.  In  fact;Cdm  nay 
decrease  to  such  a large  negative  value  that  the  overall  drag  of  the  airplane 
may  be  negative.  Of  course  to  do  this  the  discharged  air  must  be  directed 
downstream  and  be  raised  to  a momentum  level  to  sane  degree  higher  than 
when  sucked  into  the  trailing  edge.  This  last  consideration  is  not  seriously 
suggested  as  a practical  application  since  the  nower  required  to  increase 
the  momentum  of  the  discharged  air  could  no  doubt  be  more  efficiently 
utilized  in  the  airplane's  propulsive  system.  This  does,  however  exaggerated, 
indicate  some  of  the  possibilities  of  drag  control  since  the  air  must  be 
discharged  in  some  fashion,either  blown  over  seme  portion  of  the  wing  as  in 
the  Arado  System,  or  simply  thru  a discharge  duct.  If  it  is  thru  a duct  the 
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direction  could  be  varied  from  the  minimum  drag  condition  to  that  of 
maximum  drag  thereby  permitting  L/D  control  for  take  off  and  landing. 

The  drag  characteristics  of  trailing  edges  1-c  and  2-b  as  represented 
by  figures  *>1  and  52  are  erratic  in  nature.  It  is  significant  that  these 
two  configurations  represent  the  two  instances  where  a free  stream  stagnation 
point  occuired.  This,  no  doubt,  has  a great  deal  to  do  with  the  unorthodox 
appearance  of  the  drag  curves. 

Figure  82  compares  the  drag  characteristics  of  the  various  trailing 
edge  configurations.  The  drag  coefficient  is  in  this  case  the  wake  drag 
coefficient  plus  JJCq  • Trailing  edges  b-f  and  5-e  show  the 

lowest  drag  at  zero/C<^  for  the  -1°  angle  of  attack. and  h-f  also  shows  the 
lowest  drag  at  zero/C*  foi  the  9°  angle  of  attack.  Of  those  configurations 
that  reached  a jCa  - .02  practically  no  difference  in  drag  coefficient  can  be 
noted  except  .or  1-c  which  is  much  higher.  This  is  because  the  2/G^  term 
in  the  drag  coefficient  is  very  large  compared  with  the  wake  drag  except  in 
the  case  of  1-c  which  shows  a high  wake  drag  coefficient.  Of  those  tests 
that  reached  a/Cq  - .03  and  again  for /Cq  - ,0h  the  same  is  true^pointing 
out  the  great  influence  of  the  momentum  drag. 

There  is  also  an  equivalent  power  drag  which  is  a function  of  the  total 
pressure  drop  across  the  slot  and  hence  the  power  required  for  suction. 

This  will  be  discussed  in  a subsequent  section  involving  the  slot/^,or 
total  pressure  drop  coefficient 

The  effect  of  trailing  edge  suction  on  pitching  moment  is  pronounced 
and  seems  to  be  basically  a function  of  shape,  overhang  and  suction 
quantity.  The  No.  1 series  of  trailing  edges,  as  represented  by  Figures  61, 

62,  63,  and  7Qj is  characterized  by  the  pitching  moment  going  negative  with 
both  increasing  angle  and  increasing  suction.  As  the  overhang  decreases, 
there  appears  a strong  tendency  for  a discontinuity  in  the  pitching  moment 
curve  at  the  higher  angles  of  attack.  The  pitching  moments  of  the  No.  2 
and  No.  5 series  of  trailing  edges  differ  somewhat  from  those  of  the  No.  1 
series  in  that  pitching  moment  goes  negative  as  suction  is  applied  and 
positive  as  angle  of  attack  is  increased.  The  net  result  is  that  for  the 
range  ofyCQ  tested,  the  pitching  moments  are  negative  except  for  configuration 
1-c^but  the  slope  of  the  curve  is  positive.  This  is  illustrated  in  Figures 
61*,  65,  66,  and  68. 

Figure  83  is  included  to  show  a comparison  of  the  pitching  moments  versus 
for  the  several  configurations.  It  is  interesting  to  note  the  progress 
in  the  negative  direction  of  the  pitching  moments  with  the  exception  of  the 
two  trailing  edges,  1-c  and  2-b  at  13°,  which  had  the  free  stream  stagnation 
points.  For  these  two  cases  the  slope  of  the  curve  for  13-degrees  angle  of 
attack  quickly  goes  positive  after  starting  negative. 

In  general  it  appears  that  the  pitching  moment  characteristics  of 
trailing  edge  suction  are  pronounced, although  not  having  dovnwash  information 
the  true  severity  on  the  aircraft  stability  and  control  problem  cannot  be 
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estimated.  The  maximum  moment  coefficient  measured  was  -0.1*55  (Fig.  70, 
o t ■ 7.5°,>Cq  ■ 0.05, /C*  ■ 1.9).  This  value  compares  favorably  to  - 
&rr\(a.- c)"  “0»6l*  which  represents  the  pitching  moment  for  the  base  airfoil1 
equipped  with  30  percent  double  slotted  flaps  at  <T  ■ 1*0°  producing  the  same 
lift  coefficient. 

In  order  to  compare  the  various  configurations  it  is  obvious  that  lift, 
drag  and  pitch  are  but  a part  of  the  information  needed  to  gain  as  clear  an 
idea  as  possible  of  the  relative  merits  of  the  goemetry  of  the  several 
trailing  edges.  Some  measure  of  the  comparative  power  requirements  is  an 
indenspensible  part  of  such  a study.  Figures  71  through  80  show  the  drop 
in  total  press;  re  across  the  slot  as  it  is  affected  by  angle  of  attack  and 
suction  coefficient.  is  defined  as  the  total  head  drop  across  the  slot 
divided  by  the  free  stream  dynamic  pressure.  Without  suction  CCa  “ o)  the 
total  pressure  within  the  model  was  equal  to  the  free  stream  static  pressure 
and,  therefore.) different  from  the  free  stream  total  pressure  by  the  amount 
of  the  dynamic  pressure  (q).  Since,  by  definition,  this  difference  in  total 
pressure  is  divided  by  q to  obtain>Cp,  a plot  relating^/Cp  to  CX  for^C^  * o 
is  a horizontal  straight  line  (JC p ■ 1)  on  the  figures  71  to  80  regardless  of 
configuration.  Thus  Xp  is  a measure  of  the  drag  and,  therefore,  power 
requirement  of  the  slot. 

Figure  81*  is  a summary  presentation  of  the  variation  ofyCp  with^C^ 
and  trailing  edge  for  two  selected  angles  of  attack.  This  figure  shows 
that  for  a/q  » .02  and  up^  trailing  edge  1-b  suffers  the  least  total 
pressure  drop,  however,  for  lower  sCq  's  1-a  and  1-g  rival  1-b  at  least  for 
one  angle  of  attack.  The  following  tables  enable  a comparison  of  the 
efficiencies  of  the  various  slot  shapes.  Table  II  has  been  prepared  from 
Figure  81  and  shows  the  for  several  values  of  yCq  for  each  slot. 

Table  III  is  compiled  from  figure  81*  and  shows  the/Jp  values  corresponding 
to  the /\/Ojl  values  of  Table  II.  Table  IV  has  been  obtained  by  dividing 
the  values  from  Table  II  by  those  of  Table  III. 

Table  II 


at  13°  angle  of  attack  for  several  values  of 


1-a 

1-b 

1-c 

1-g 

2-b 

2-d 

2-e  l*-f  5-e  -a 

.01 

.61* 

.63 

.06 

.36 

.08 

.76 

.1*1  .08  .12  .22 

.02 

1.28 

1.05 

co 

o 

• 

CM 

r- 

• 

-.05 

1.19 

.61  1.07  »2li 

.03 

.30 

1.73 

.89 

.01* 

1.1*6 

1.81* 

.99 
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Table  III 


/Cp 

at  13  degrees 

angle  of 

attack  for  several  values  of  /Cq 

A? 

l-a 

1-b 

1-c 

1-g 

2-b 

2-d 

2-e  l*-f 

5-e  d-a 

.01 

1.9 

l.l* 

1.5 

1.1* 

3.1* 

2.8 

2.7  3.1* 

5.0  1*.0 

.02 

3.6 

2.0 

2.5 

2.7 

10.0 

8.0 

.03 

6.9 

3.3 

6.7 

lu5 

.04 

5.1* 

10.5 

7.0 

Table  IV 


at  13 

degrees 

angle 

of  attack 

for 

several 

values 

of/$* 

c 

l-a 

1-b 

1-c 

1-g 

2-b 

2-d 

2-e 

l*-f 

5-e  d-a 

.01 

.336 

.1*50 

.01*0 

.258 

.023 

.272 

.152 

.021* 

.021*  .055 

.02 

.358 

.525 

.032 

.267 

.119 

.131* 

.03 

.01*3 

.525 

.198 

.01* 

.31*1 

.11*1 

Table  IV  indicates  that  trailing  edge  configuration  1-b  was  the  most 
efficient  one  tested  with  respect  to  power  required  to  move  the  air  through 
the  slot.  A more  significant  comparison  of  the  efficiencies  is  presented 
in  figure  81*-a  in  which  suction  coefficients  versus  total  drag  coefficients 
are  plotted  for.  an  angle  of  attack  of  90°.  The  total  drag  coefficient 
has  been  determined  by  the  addition  of  the  profile  drag  (/Cdn) ) to  the 
equivalent  power  drag  coefficient  (see  section  IV  B).  Such  a definition 
involves  two  assumptions  i 

1.  That  the  momentum  drag  be  zero. 

2.  That  the  suction  system  except  for  the  slot  have  an  efficiency  of 
100  percent.  ' 

The  profile  drag  coefficient  of  the  base  airfoil  (unmodified  NACA  23015) 
at  9-degrees  angle  of  attack  was  determined  to  be  0.02.  This  value  appears 
as  a horizontal  straight  line  on  figure  81*-a  for  purposes  of  comparison.  It 


is  evident  from  this  plot  that  configuration  1-b  is  again  superior  since  it 
shows  the  minimum  drag  coefficient  (0.011)  and  also  since  it  shows  the  greatest 
spread  of/63  (.005  to  .018)  which  yields  a total  drag  coefficient  less  than 
that  of  the  base  profile  at  this  angle. 

The  Variable  Overhang  Trailing  Edge 

Of  the  several  geometric  variables  common  to  the  types  of  trailing  edge 
suction  slots  studied  the  amount  of  lower  lip  overhang  seems  to  be  one  of 
the  most  significant.  For  this  reason  a variable  lower  lip  was  constructed 
and  tested  in  combination  with  upper  lip  number  one  (see  Figure  33).  Figure 
85  shows  a comparison  of  the  lift  coefficients  obtained  vs. yC^  for  each  of 
several  angles  of  attack  of  the  different  amounts  of  overhang  tested.  The 
notations  1-v-l  to  l-v-10  designates  increasing  overhang  of  the  lower  lip. 
Several  interesting  phenomena  are  apparent  from  these  plots: 

1.  As  the  lower  lip  is  extended  there  is  a gradual  decrease  in  the  lift 

coefficient  for  the  angle  of  11  degrees  at  every  value  of /Cq  including jCq  » o . 

2.  As  the  lower  lip  is  extended  there  is  a decrease  in  lift  coefficient 

for  the  angles  -1  and  5 degrees  for  most  all  values  of  jcq  with  the  exception 
of/c^  ■ o. 

3.  As  the  lower  lip  is  extended  the  dip  which  occurs  in  the  -1°  lift 
curve  tends  to  flatten  and  this  curve  finally  appears  as  a straight  line. 

Figure  86  shows  an  idealized  explanation  of  the  boundary  layer  removal 
at  a low  angle  of  attack  and  with  a relatively  short  overhang  of  the  lower 
lip.  Without  suction  at  01  -0°  or  -1°  there  is  a significant  boundary  layer 
thickness  on  both  surfaces  of  the  profile.  As  suction  is  added  the  upper 
boundary  layer  tends  to  be  removed  first^having  sn  effect  upon  the  lift 
comparable  to  an  increase  in  camber,  thus  producing  a lift  coefficient 
corresponding  toyc*  « A.  As  the  suction  is  increased  to  B the  trailing 
edge  sink  apparently  has  less  effect  upon  the  free  streamlines  than  upon 
the  boundary  layer  on  the  lower  surface.  By  removing  a large  portion  of  the 
lower  surface  boundary  layer  the  effect  upon  lift  is  comparable  to  a decrease 
in  camber  near  the  trailing  edge.  This  results  in  the  dip  of  the/C vs. 

curve  having  a minimum  at/C $ • B.  Beyond  this  suction  quantity  the 
free  streamlines  above  the  trailing  edge  of  the  profile  are  more  and  more 
affected  with  increasing  suction.  As  the  angle  of  attack  increases  one  would 
expect  a thickening  of  the  upper  surface  boundary  layer  and  a reduction  in 
thickness  of  the  lower  surface  boundary  layer.  This  would  explain  the 
absence  of  such  a dip  in  the  curve  st  higher  angles.  Also  as  the  lower  lip  is 
extended  there  would  be  less  tendency  for  flow  around  the  trailing  edge  into 
the  slot, which  would  explain  the  lack  of  the  dip  in  the  curve  for  the  greater 
overhangs.  This  discussion  referring  to  figure  86  is  offered  as  a suggested 
explanation  for  the  dip  in  the>C^  versus  curve  as  shown  in  Figure  85, 
trailing  edge  1-v-l,  for  an  angle  of  attack  of  -1  degree.  This  shspe  of 
curve  has  been  observed  in  other  tests  made  with  the  force  model  and  pressure 
model,  the  detailed  data  of  which  are  presented  in  reference  6. 


hi. 


Figure  87  shows  the  AO  , for  several  values  of  •nd  angle  of 
attack,  plotted  against  the  ratio  a/ cy  where  a is  the  chord  of  the  upper 
surface  and  c is  the  chord  of  the  lower  surface  of  the  force  model  fitted 
with  the  variable  trailing  edge.  This  illustrates  the  rapid  increase  of 
lift  coefficient  at  a given  value  of /Cq  with  increasing  values  of  a/c. 

This  of  course  agrees  with  theory  which  calls  for  an  increasing  as  the 

slot  approaches  the  trailing  edge.  It  is  interesting  to  note  that  at  the 
highest  angle  shown  (10  degrees)  there  is  an  optimum  value  for  this  ratio  a/c 
for  two  values  of  and  evidently  this  optimum  value  increase  with  • 

Comparison  of  Results  with  those  of  other  Experimenters 

Figures  88  and  89  coir, pare  some  of  the  results  herein  reported  with  these 
reported  by  several  other  experiments lists.  The  comparison  made  in  Figure  88 
is  with  the  optimum  trailing  edge  slot  location  reported  by  United  Aircraft1 2 3. 
Trailing  edge  configuration  5-e  was  chosen  for  this  because  it  produced  the 
greatest  A Cjl  at  zero  angle  of  attack.  It  appears  to  be  somewhat  more 
effective  than  the  U.A.C.  wing,  at  this  /Ca  (.01),  but  this  may  be  due  to 
the  difference  in  sections.  Figure  89  compares  trailing  edge  1-b  with  those 
of  Ehlers*  and  SegenscheiV.  All  three  represent  tests  made  with  the  same 
basic  section  (N.A.C.A.  23015)  but  with  somewhat  different  trailing  edge  slots. 
Configuration  1-b  seems  to  be  inferior  at  the  lower  angles  of  attack  but  from 
ot  ■ 7.5  degrees  and  up  it  is  more  effective  than  those  reported  by  both 
Ehlers  and  Regenscheit.  ^A'Coi  » .01  has  been  selected  for  this  comparison 
because  it  is  felt  that  it  is  in  the  range  of  practical  application  and 
interest. 

Conclusions 

It  seems  possible  to  deduce  the  general  behavior  of  subsonic  profiles 
equipped  with  trailing  edge  suction  slots  fran  this  study.  Hie  characteristics 
of  lift,  drag,  pitching  moment,  power  requirement,  slot  shape  and  general 
geometry.as  influenced  by  suctionyhave  been  investigated  in  sane  detail  and 
it  is  felt  that  the  following  statements  are  justified  for  sn  airfoil  without 
a flap. 

1.  Hie  lift  coefficient  of  a profile  equipped  with  a trailing  edge 

suction  slot  can  be  increased  by  the  application  cf  suction  except 

for  the  case  where  a free  stream  stagnation  point  occurs;  in  such  an  instance 
the  lift  increment  can  become  negative. 

2.  Hie  lift  coefficient  of  such  an  airfoil  in  general  increases  as  the 
sink  approaches  the  trailing  edge  as  shown  by  theory  and  experiment.  Hie 
exception  to  this  statement  is  the  optimum  value  of  the  ratio  a/c  beyond 
which  the  opposite  is  true. 


1. ~  ’TSTerence"  IE)' 

2.  Reference  11  Fig.  19 , Page  39 

3.  Reference  12  Fig.  12,  Page  29 


3»  The  lift  coefficient  is  increased  at  first  (low  values  of ) 
by  boundary  layer  removal  which  permits  a flow  more  nearly  approaching  that 
indicated  by  potential  flow  theory.  At  higher  values  of  /Oq  the  lift  increase 
can  be  thought  of  as  a result  of  a circulation  increase  and  sink  effect. 

b.  The  lift  increase  with/C*^  is  very  sensitive  to  the  general 
geometry  cf  the  trailing  edge  slot  but  it  dees  not  seem  at  all  sensitive  to 
the  shape  of  the  upper  lip  inlet  contour  unless  this  shape  is  extreme. 

5.  The  drag  of  such  a configuration  is  variable  depending  upon  how 
the  removed  air  is  returned  to  the  free  stream.  The  wake  may  be  practically 
eliminated  with  the  application  of  enough  suction. 

6.  Pitching  moment  is  in  general  negative  and  moderate  increasing  in 
the  negative  direction  with  a positive  increase  of  /Cq  . 

7.  The  power  requirement  can  be  represented  by  the  total  pressure  loss 
across  the  slot  and  it  is  a function  chiefly  of  the  slot  width  which  determines 
slot  velocity  for  a given  /Cq  . It  is  also  a function  of  angle  of  attack 

and  slot  lower  lip  overhang. 
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'DEAliZED  EXPLANATION  OF  BOUND A R < LAVER  REMOVAL 
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DARY  LAVER  REMOVED  FROM  UPPER  SURFACE 
U"f  STIUL  E/lSTlNS  CM  LOWER  SURFACE. 
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